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Abstract 

We develop techniques for studying the effects of self-interactions in the conformal 
sector of an unparticle model. Their physics is encoded in the higher n-point functions 
of the conformal theory. We study inclusive processes and argue that the inclusive 
production of unparticle stuff in standard model processes due to the unparticle self- 
interactions can be decomposed using the conformal partial wave expansion and its 
generalizations into a sum over contributions from the production of various kinds 
of unparticle stuff, corresponding to different primary conformal operators. Such 
processes typically involve the production of unparticle stuff associated with oper- 
ators other than those to which the standard model couples directly. Thus just as 
interactions between particles allow scattering processes to produce new particles in 
the final state, so unparticle self-interactions cause the production of various kinds 
of unparticle stuff. We discuss both inclusive and exclusive methods for computing 
these processes. The resulting picture, we believe, is a step towards understanding 
what unparticle stuff "looks like" because it is quite analogous to way we describe 
the production and scattering of ordinary particles in quantum field theory, with the 
primary conformal operators playing the role of particles and the coefficients in the 
conformal partial wave expansion (and its generalization to include more fields) play- 
ing the role of amplitudes. We exemplify our methods in the 2D toy model that we 
discussed previously in which the Banks-Zaks theory is exactly solvable. 
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1 Introduction 



Since the original formulation of unparticle physics [HE] as an effective field theory in which 
the standard model couples only at high energies to a "Banks-Zaks" [3j HJ E] sector that 
is scale-invariant at low energies, this idea has been explored both theoretically and phe- 
nomenologically in many papers. We start by mentioning what we find the most interesting 
theoretical developments. 

The way unparticle physics arises from weakly coupled Banks-Zaks- like theories was 
analyzed explicitly in several examples in [6j |7J. Strongly coupled Banks-Zaks sectors can 
include theories in the conformal window of supersymmetric QCD [8] as discussed in [91110116]. 
or other supersymmetric or non-supersymmetric gauge theories [TT1 IT2l [T3] . Quite generally 
(no counterexamples are known), unitary scale- invariant theories have the full conformal 
invariance [T41 [15]. This imposes lower bounds on operator dimensions and dictates the 
tensor structure of the unparticle propagators [7J [16j [17]. The dimension of a primary^ 
vector operator (9 M must be as large as du > 3. This would usually suppress the possible 
interactions with the standard model both in the absolute magnitude and relative to the 
accompanying standard model contact terms |7j. Similarly, for an antisymmetric tensor 
0^ v du > 2 and for a symmetric traceless tensor 0^ u du > 4. However, scalar operators 
can start from d u > 1. Spin-| operators must have du > 3/2, and they mostly become 
relevant if it is possible to describe unparticle stuff charged under the standard model gauge 
interactions (otherwise, the fermionic unparticle operator can only couple to the standard 
model fields similarly to a singlet neutrino.) Conformal invariance also introduces a specific 
(^-dependence m t° the tensor structure of the 2-point functions in momentum space [71 [T7]. 
It is also interesting to mention the work [18] that used the conformal symmetry to study 
energy-momentum correlations in unparticle stuff. In particular, they find that the total 
momentum flux in a particular direction is always proportional to the energy flux, like for a 
massless particle, and derive the angle-dependent correlation functions of energy and charge. 

The possibility that the coupling of an unparticle operator to the Higgs would break the 
low energy conformal symmetry of the Banks-Zaks sector was discussed in p]. An interac- 
tion of the form (l/M^ uv ~ 2 )\H\ 2 O uv that flows in the IR to (A^" d "/M^ v " 2 )|i7| 2 C w 
becomes a relevant operator in the conformal sector (for du < 4) introducing a scale 
Aj" d " ~ (Au/M u ) duv ~ du M u ~ du v 2 , where v is the Hiees vev. This can break the conformal 
symmetry at energies E < Ay (while preserving unparticle behavior for <C E <C A^). 
The consequences of the breaking depend on the particular realization of the Banks-Zaks 
sector, and in many cases rich and surprising phenomenology is expected (see [T9l 1201 l2Tj 
and especially the hidden valley picture of [6]). 

The fact that some conformal field theories are dual to gravitational theories in anti-de 
Sitter space with one extra dimension (the AdS/CFT correspondence [221 E3J EI] or its more 

*A primary operator is one that behaves covariantly under conformal transformations. Any operator of a 
definite scaling dimension that is not a derivative of another operator is primary. In 2D CFTs these operators 
are often referred to as "quasi-primary," while the term "primary" refers to the properties of operators with 
respect to the full Virasoro algebra. 
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phenomenological variants, such as the Randall-Sundrum models [251 [26]) can be useful for 
studying various aspects of unparticle physics from a different perspective. For example, 
some works on unparticle phenomenology involving internal scalar unparticle propagators 
obtained divergent results for d u > 2 [271 1211 [29] . Looking from the AdS perspective, the 
authors of [30] (following an earlier work |31j) explained that for du > 2 (or du > 5/2 for 
fermions) the propagator includes UV-dependent terms that do not decouple when the UV 
cutoff is removed. These terms cancel the divergences, and in fact dominate much of the 
physics in this range of du- As pointed out already in [3 lj (see also [32]), these counterterms 
are needed in order for the two-point function to be a well-defined distribution. In the 
context of unparticle physics, these are just the standard model contact terms generated 
from integrating out the high energy physics that couples the standard model to the Banks- 
Zaks sector. An explicit example of this was presented in [7J. The renormalization group 
running of such terms has been analyzed from the holographic perspective in [33]. On the 
other hand, the UV completion provided by the RS2 model [25J has been analyzed in [34J. 
The AdS/CFT correspondence was also used in [30] for analyzing how the unparticle stuff 
can couple to gauge interactions, similarly to their earlier suggestion of gauging a non-local 
action [35] (see also [361 S3 EH [39]). Other uses of the AdS picture included describing 
Banks-Zaks sectors in which the conformal invariance is broken in the IR [lOl El [30] and 
analyzing the possibility suggested in [H] that "unhiggs" is responsible for the electroweak 
symmetry breaking [4"2] . 

The disadvantage of the AdS/CFT-based approaches is the difficulty to explicitly specify 
the 4D Lagrangian description of the physics. The original AdS / CFT setup can describe the 
Banks-Zaks sector in its low-energy limit, the CFT. However, since the number of examples 
from string theory is limited, typically one needs to pick the field content and the Lagrangian 
in the AdS space "by hand" p without knowing what 4D theory is being described. Many 
CFTs do not have weakly coupled AdS duals at all. Those that do are typically large-iV 
gauge theories with large 't Hooft coupling A (see also [13]). There is no special reason to 
believe that the conformal sectors in our world should belong to this class of theories. Some 
properties may work out anyway, as in AdS/QCD [S1H5], but it can often be hard to separate 
a real effect from an artifact of an uncontrolled approximation. Certain features of unparticle 
physics that are missed by assuming N — > oo, A — > oo were discussed in [6]. Furthermore, 
while knowing correlation functions of CFTs is important for unparticle physics, one would 
like to include the coupling to the standard model as well and describe the UV completion 
of the combined theory. In order to add these non-conformal ingredients, the AdS space 
needs to be modified in the IR. For example, the RS2 model [25] cuts off the AdS space 
by a brane. As a result, the CFT gets cut off in the UV due to interactions with the 
boundary values of the AdS fields (gravity and others) which become dynamical fields in the 
4D theory [HU H7J [31]. Unparticle physics aspects of an RS2 model with a massive vector 
field in the bulk were analyzed in [31] . In general, different ways of modifying the AdS space 
in the IR and regularizing its physics will correspond to different choices of the coupling 

2 It is widely believed that knowing how (or whether) such a choice can actually be realized in string 
theory is not essential in order for the correspondence to work. 
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to the standard model in the UV. But unfortunately the relation between the two sides of 
this extended correspondence is far from being straightforward. Therefore, while AdS-based 
models can provide very useful guidance and examples, their ability to describe realistic 
unparticle physics scenarios is limited. 

The understanding of unparticle physics is incomplete without taking into account the 
self-interactions of the low energy conformal sector. The main goal of the present paper 
is to contribute to this understanding. So far most works have focused on the 2-point 
function of the unparticle operator O to which the standard model couples, with a only 
few excursions [6j HE] into the more complicated higher n-point functions which contain the 
information about the interactions. 

Scale invariance requires the 2-point function of O of dimension d to have the forrrJ§ 

(0|TO(x)O(0)|0> ex - 1 , oc / (-p 2 - ie) d ~ D/2 (1.1) 

where D is the spacetime dimension. There are two ways in which this 2-point function can 
appear in physical processes. First, it can appear as an internal line in a process that includes 
two standard model-unparticle interaction vertices. The momentum-space expression for the 
2-point function makes the calculation straightforward. This results in interesting effects due 
to the non-zero imaginary part at all p 2 [21 127] . The other physical effect is the production 
of unparticle stuff. The corresponding phase space can be determined either from scale 
invariance [T] or by computing the imaginary part of the 2-point function [2]: 

$oc (p 2 ) d ~ D/2 6(p°)6(p 2 ) (1.2) 

These two types of processes describe the physics at the leading order in the standard model- 
unparticle coupling. 

Higher order processes require additional machinery since they will depend on 3- and 
higher n-point functions of O. Furthermore, we will argue that it is useful to not restrict 
our attention solely to the operator O that couples to the standard model sector. Instead 
we will use more of the power of conformal invariance and consider the primary operators 
Oj of dimension dj of the conformal field theory. Conformal invariance requires the 2-point 
functions to have the form@ 

(0|T^.(x)O fc (0)|0) cx 5jk ex S jk ( T^e"^ {-p 2 - ie) d ~ D ' 2 (1.3) 

(— x l + le) J 

The Sjk in f ll .3D allows us to identify different primary operators with different kinds of 
unparticle stuff. For each Oj, there is a unique phase space given by 

$j<x(p 2 ) dj - D/2 9(p°)9(p 2 ) (1.4) 



3 For simplicity of presentation, we assumed here that the operator is a Lorentz scalar. 

4 If there is more than one operator with a given dimension and tensor structure, we can choose linear 
combinations to get the 5jk- Again, for simplicity of presentation, we assumed here that the operators 
are Lorentz scalars, but most of what we say can be easily generalized to higher tensors. Various general 
properties of conformal theories are discussed in [321 EH [SH E21 E] • 
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The form of the 3-point function is also fixed by conformal invariance. In particular (for 
example), for primary scalars Oi with dimensions df. 

(0\TO 1 (x 1 )O 2 (x2)O 3 (x 3 )\0) oc , . ^ . \ —: wo . T—r—; wo (1.5) 

\ X 12J \ X 13J l x 23/ 

where OC ij — OC >i OC j . This was used in [48J (in the case where all Oi are the same) to 
study processes that involve 3 internal standard model-unparticle vertices. Higher n-point 
functions are highly constrained by conformal invariance, but not completely determined. 

Our goal in this paper is to understand how to use this structure to analyze processes in 
which unparticle stuff is produced as an outgoing state. We will discuss two approaches. We 
will show how to analyze inclusive processes by looking at the discontinuities across physical 
cuts in the n-functions of O. And we will argue that these discontinuities are related to a 
sum over primary conformal operators of squared "amplitudes" for the production of the 
corresponding unparticle stuff. These amplitudes, in turn, are determined by the coefficient 
functions in the conformal partial wave expansion. The resulting picture, we believe, is a 
step towards understanding what unparticle stuff "looks like" because it is quite analogous to 
way we describe the production and scattering of ordinary particles in quantum field theory, 
with the primary conformal operators playing the role of particles and the coefficients in the 
conformal partial wave expansion (and its generalization to include more fields) playing the 
role of amplitudes. 

We will present the basic ideas, which are valid for any conformal theory in any number 
of dimensions, in section |2j and dedicate most of the rest of the paper to testing them on the 
2D example of unparticle physics that we discussed in [53]. The Banks-Zaks sector in that 
example is the Sommerfield model of massless fermions coupled to a massive vector field. 
This model is exactly solvable and flows to the Thirring model in the infrared, as we will 
discuss in section |3l We couple the Sommerfield model to a toy standard model, which is 
simply a massive scalar carrying a global U(l) charge. In the infrared, the resulting inter- 
action flows to a coupling of two charged scalars to an unparticle operator with a fractional 
anomalous dimension. In section H] we apply the operator product expansion to the solution 
of the Sommerfield model to find the exact 2n-point functions of the unparticle operator. 
We discuss the mathematical structures that appear in the solution and certain infrared 
issues involved in looking at them in momentum space. Before studying self-interactions, 
we briefly review in section [5] the simplest unparticle process shown in figure [T^l in which 
two toy standard model scalars "disappear" into unparticle stuff. Because we have the exact 
solution for the Banks-Zaks sector correlation functions, we can see precisely how the system 
makes the transition from the low-energy unparticle physics to the high-energy physics of 
free particles. The answer, as we discussed in [53], is rather simple. The "spectrum" of the 
model that we can see in that process consists of unparticle stuff and massive bosons. As the 
incoming energy of the standard model particles is increased, the unparticle stuff is always 
there but more and more massive bosons are emitted and the combination becomes more 
and more like the free-fermion production cross-section. In sections [6] and [7] we analyze the 
missing charge and missing energy processes shown in figures [lb and [It, respectively. These 
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Figure 1: (a) A disappearance process, (b) A missing charge process, (c) A missing energy process. 
The arrows indicate the flow of the global U(l) charge. 



inclusive processes are mediated by the conformal sector self-interactions. We show both 
inclusive and exclusive methods of calculation. The spectrum again includes unparticle stuff 
and massive bosons, but the unparticle stuff corresponds to operators other than those to 
which the standard model couples directly (there are also massless excitations that corre- 
spond to operators of integer dimensions). We comment on the massive bosons in section |HJ 
In section [9] we summarize the conclusions. 
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2 Unparticle self-interactions 



A useful concept in conformal field theories is the conformal partial-wave expansion [511 155| 
[561 ETJ EH [52] (see also [58j[59j[60] and references therein in the context of AdS/CFT), which 
is a generalization of the operator product expansion (OPE). For any two operators Oi(xi) 
and 2 (x 2 ) with an arbitrary separation between them, it is possible to write 

TO 1 (x 1 )O 2 (x 2 )\0) = Ysf d D xiQ k (x\x 1 ,x 2 )O k (x)\0) (2.1) 

k 

where T denotes time-ordering, O k are the various primary operators in the theory (unlike in 
the OPE, their derivatives need not be included separately) and the coefficients iQ k (x\xi, x 2 ) 
are the 3-point functions of Oi(xi), 2 (x 2 ) and O k (x), with the Ok leg amputated, namely 

J d D x (0\TO k (x')O k (x)\0}tQ k (x\x u x 2 ) = (0\TO k (x')O 1 (x 1 )O 2 (x 2 )\0) (2.2) 



If 0\ and 2 are scalars, the operators Ok are completely symmetric traceless tensors [57] 
whose amputated 3-point functions Qk{%\x\,x<i) are fully determined, up to a constant pref- 
actor, by the dimensions of the three operators and the tensor rank of Ok- 

Amputated 3-point functions are exactly what is needed for computing processes that 
produce unparticle stuff corresponding to the operator whose leg is amputated. This follows 
because using (12. ip twice we can write 

(0\TO* 2 (x 2 )O* 1 (x 1 )O 1 (y 1 )O 2 (y 2 )\0) 

= J2[ d D xd D yQl(x\x 1 ,x 2 )(0\TOl(x)O k (y)\0)Qk(y\yi,y2) (2.3) 

For example, suppose we had the coupling 

Ant oc <j) 2 (2.4) 

where is a standard model field and O is an unparticle operator. According to ( 12.3H with 
0\ = 2 = O, by taking the discontinuity across the cut of the 4-point function of O 
(figure [2]), which corresponds to the inclusive process 

+ — > <fi + <f) + unparticle stuff (2.5) 

we obtain the sum of cross-section for the processes 

+ 0^0 + + {O k stuff} (2.6) 

where Ok are the various primary operators in the theory (that do not necessarily couple 
to the standard model directly). The amplitudes M. of these processes are the amputated 
3-point functions Qk (times factors coming from the standard model), while cutting the Ok 
propagators in (12. 3p gives the O k phase spaces in [\1A\ . 



Figure 2: On the left: the unparticle 4-point function (schematically) with a cut over the unparticle 
stuff for computing the cross-section of the process (12.51) , inclusive with respect to the unparticle stuff. 
The dashed lines are the standard model particles 4>. Each vertex is the interaction (j2.4p . (One should 
include also diagrams in which the standard model particles are attached to the 4-point function in other 
possible ways, as we will do in section 0) On the right: representation (12. 3D of the 4-point function as a 
sum of terms, each with two amputated 3-point functions Qj, connected by an Ok propagator. The cut 
through the Ok propagator allows us to interpret the terms in the sum as related to the cross-sections 
of the exclusive processes (j2.6l) . 




Figure 3: The D'EPP formula in its original form (12.71) (top) and in the form (12.81) that is useful for 
amputation (bottom). 

There are powerful tools in conformal field theory for working out these amplitudes 
explicitly. The amputated 3-point function is formally proportional to the ordinary 3-point 
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function of 0\ and 2 with a "shadow operator" corresponding to the third operator Ok- 
The shadow operator has the same tensor rank i as Ok, but dimension A = D — A, where 
A is the dimension of Ok and D is the spacetime dimension. This method is based on the 
D'EPP formula [61] (see also [Ml EH 1521 E2]), which says that for S 1 + 5 2 + 5 3 = D (in 
Euclidean space) 



f d n x rgp r{6 2 ) t(6 3 ) = nD/2 r(f-^)r(£-^r(f-jO 

/ / 2 / 2 \fe / 2 / 2 \D/2—&3 1 2 \D/2—&2 I 2 \D/2—S\ 

J \ x u) \ X 24J \ X 43) \ X 12) \ X 13) \ X 23) 

For solving (12.21) we can use this formula as 

1 fjn_ [r(A/2)] 2 r(A) D/2 T(D/2 - A) [T(A/2)] : 



(2.7) 



" X4 . „ . A /o . „ . A /o . „ . A " , n xrf-A/2 / n s A/2 / n s A/2 



(* 2 2 ) d ~ A/2 ■> (x? 4 ) A/2 (x| 4 ) A/2 (x| 3 ) A ' (x? 2 ) d ~ A/2 (x? 3 ) A/2 (^ 3 ) A/2 

This is described graphically in figure El Since the 3-point function (II. 5ft of two scalars 0\ 
and 2 of dimension d with a third scalai0 Ok of dimension A, and the two-point function 
of Ok, are 

WO^)0 M{x3m = ( ^ 2 + g) ^ ( _^ ;e) , /2( _ x|3 + g) . /2 («> 

<0|TO t (z)O»(0)|0) = . . (2.10) 

(— X z + 16) 

where C3 and C 2 are constants, we obtain the amputated 3-point function to be 

[r(A/2)] 2 r(A)c 3 

n D/2 [r(A/2)] 2 T(D/2 - A) C 2 

X (-x? 2+ ^- A/2 (-x? 4 + ,e) A/2 (-x| 4 + , e ) A / 2 (2 ' 11} 



Q fe (x 4 |xi,x 2 ) 



In momentum space it gives the amplitude 

2 D - M+A r(A)r(D-rf- a/2)c 3 , . . 

M - T { D/2-l)r(d + A,2-D/U I{P ' Q) <2 ' 12) 
(times the standard model factors), where 

/dPk 

(-P - % ef- d -^ HP - ^ - .e) A / 2 (-(A: - (P - Q)) 2 - *e) A/2 (2 ' 13) 

where P is the momentum incoming from the standard model at x\ and Q is the unparticle 
momentum outgoing at x 4 . The phase space of Ok is 



7T D / 2+1 C 2 



2 2A-jj-i p(A) T(A - D/2 + 1) 



(Q 2 ) A ~ D/2 fl(QV(Q 2 ) (2.i4) 



5 A generalization for higher tensors exists as well. We will consider a vector operator in section IT^fl 
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Note that $ is positive if C2 > and A > D/2 — 1, which is exactly the well-known bound 
on the dimension of a scalar operator [TTl [52] . For the limiting case A = D/2 — 1, use 



with e = A - (D/2 - 1) to obtain 



o D/2+1 (-1 

HQ) = T(D/2-l) 6iQ0)6{Q2) (2 - 16) 



We thus have 



■ M{ 2 ~ _ ^ D - Ad+1 * D,2+l r(A) [T(D - d - A/2)] 2 C 2 

1 1 [T(D/2- A)] 2 [T(d + A/2- D/2)} 2 T(A- D/2 + 1) C 2 

x\I(P,Q)\ 2 (Q 2 ) A ~ D/2 0(Q°)9(Q 2 ) (2.17) 
We can combine the denominators in (12.131) by using 



n r K' 

j 



which gives 



9$ ^ 2 3 ^+ 1 ^/ 2 + 1 r(A)[r(D-rf + A/2)] 2 a 2 



[T(D/2 - A)] 2 [r(d + A/2 - D/2)] 2 [r(A/2)] 4 T(A - D/2 + 1) C 2 

x 



I(P,Q) 2 (Q 2 ) A - D/2 9(Q°)9(Q 2 ) (2.19) 



where 



/(P, Q) = dat\ da.2 da 3 5 j a, — 1 J 



D-d-A/2-1 A/2-1 A/2-1 



a 2 a 3 



2.20) 



(P + ,(P,Q) + , e ) D - d+A / 2 
where we shifted the integration variable k and defined 

g(P, Q) = a x {l- a^P 2 - 2a x a z QP + a 3 (l - a 3 )Q 2 (2.21) 

In the following sections we will do the inclusive and exclusive calculations explicitly in 
a 2D toy model, the Sommerfield model, where we consider the processes whose inclusive 
descriptions are given in figures [Tb and c. 
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3 2D toy model of unparticle physics 



3.1 Sommerfield model of a Banks-Zaks sector 

In this paper we apply the ideas of section 2 to analyze the unparticle stuff that appears 
in the Sommerfield model [63], EU |65] 166 ] [67] . that is the Schwinger model [68] with an 
additional mass term for the vector bosonO 

1 m 2 

C = iP{ift -e4)iP--F^F^ + ^A^ (3.1) 

We are interested in this theory since, like the Schwinger model, it is an exactly solvable 
model that becomes scale-invariant at low energies, and (unlike the Schwinger model) has 
fractional anomalous dimensions. At low energies, below the vector boson mass, the theory 
reduces to the Thirring model, that is the theory of a fermion with a quartic self-interaction 
which is scale- invariant at all energies [691 Ell EE]- 

In order to solve the model, it is convenient to decompose sJ^ 

A» = d»V + e^d v A (3.2) 

The Lagrangian becomes 

1 777^ 

C = iip flip- eifopil) (d»V + e^duA) + -A a 2 A + {d^Vd^V - d^Ad^A) (3.3) 
If we change the fermionic variable to 

^ = e ie ( v+ ^ 5 )v (3.4) 

the fermion becomes free: 

2 i 2 
777 777 

C = m + -2^V^V + -A a 2 A - —d^Ad^A (3.5) 
In the last term of (I3.5P ml has been replaced by 

e 2 

m 2 = ml H (3.6) 

7T 



6 Our conventions, as in [53], are: g 00 = —g 11 = 1, e° = — e 10 = — eoi = eio = 1- From the defining 
properties {7^,7"} = 2g» v and 7 5 = -\e^Yl w , it follows that j^ 5 = -e^7„ and 7^7" = g pv + e^7 5 , 

and we will use the representation 7 = ^ J , 7 1 = f ^ ^ , 7 s = 7°7 1 = ( J ^1 ) Then the 

components ip\ and "02 describe a right-moving and left-moving fermion, respectively. 
7 In 2D, an arbitrary vector can be expanded in terms of any non-null vector fc M as 

_ k^kP - e^K^K 
~ P p 

In the Schwinger model (mo = 0), V would be the unphysical longitudinal polarization that can be set to 
zero in the Lorenz gauge. 
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in order to account for the fact that the path integral measure is not invariant under the A 
part of Q23) [72] @ 

We can now calculate fermionic n-point functions as 

(0\Tip a (x) ...ip;(y)...\0) = (0|Te^ e ( v(:E)+ ' 4{:c)7 ') . . . e ie ( v(y)+ - 4( ^) . . . |0) 

x{0\T* a (x)...*}(y)...\0) (3.7) 

where the subscript in 7^ specifies that it acts on the spinor at point x (in our representation, 
7 5, 0i = 7 5 V'2 = — ip2-) Using Wick's theorem and the observation that 



-% j d 2 xe ipx (0|TV(a;)V(0)|0> 
-i [ d 2 xe ipx (0\TA(x)A(0)\0) 



m\ p 2 



1 



2\2 



m 2p2 



\p 



(3- 

J2 ) (3-9) 



we obtain that the n-point function is given by the corresponding free n-point function 
multiplied by 



(3.10) 



j>i 

where i, j run over the n points and 

^2 



C (x) 
C{x) 



exp 
exp 
exp 



rat 



[D(x) - D 



oc 



(-x 2 + ie) 



— e 2 /47rm,Q 



(3-11) 



i^-[(A(x)-A(0))-(D(x)-D(0))} 



2iTm 2 - 



K [mV—x 2 + ie ) + In ( (mV-i 2 + 



ie 



with 



where we useql 



olE 



A(x) 
D(x) 



d 2 p 



-ipx 



(2ir) 2 p 2 — m 2 + ie 
d 2 p e~ ipx i 
(2ir) 2 p 2 + ie An 



■■^-Ko [rn\/—x 2 + iej 



In 



— x + ie 



(3.12) 
(3.13) 

(3.14) 
(3.15) 



8 The same effect gives mass e/y/ir to the gauge boson in the Schwinger modei. See also [73] . 
9 Kq is the modified Bessel function of the second kind and xq is an arbitrary constant that will cancel 
out in the following. For y — > 00, K (y) ~ rE. e -v — > 0. For y — > 0, K a (y) — — ln(y/2) — + 0(y 2 ), where 

7B = -T'(l) ~ 0.577 is Euler's constant. Note that A(0) - £>(0) = (i/2n) In (e 7E x m/2) = (i/2n) In (£x m) 
in C{x) is finite. The D(0) term in Cq(x) does not have any effect since any n-point function has —n/2 such 
factors so their presence is equivalent to changing the normalization of the fermion field. 
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and r/ij and Kij are sign factors that depend on whether the operators at the two points are 
the fields or their adjoints and their handedness: 

n ^and^* K ..=^5_f+1, « = P ( o m) 

VlJ \ -1 , (</> and V) or (</>* and J ~ 7l 7j \ -1 , ff^ 1 J 

Similar expressions for the correlation functions have been obtained by various methods in 
the past [31IZSICS]. 

In the short- distance limit, (xi — Xj) 2 <C 1/m 2 , C(x) — > 1 and one obtains free-fermion 
behavior q In the large-distance limit, (xi — Xj) 2 1/m 2 , K does not contribute, and C{x) 
is just a power of x 2 \ 

C( x )^ [(( ffl ) 2 (-, 2 + if 2/W (3.17) 

leading to scale-invariant behavior with fractional anomalous dimensions. This is the unpar- 
ticle regime. Note that m plays the role of i\u from pQ. 

Additional properties of the Sommerfield model are described in appendix |A] 

3.2 Coupling to the "standard model" 

We assume that the very high energy theory includes the interaction 

Ant = | ^(1 + 7 5 )X0*+?(1-7 5 )X0] +h.c. 

= ^(V 2 *Xi0* + ViX2 0)+h.c. (3.18) 

that couples the fermion ip of the Banks-Zaks sector to a neutral complex scalar with mass 
<C m that plays the role of a standard model field. The interaction is mediated by the 

heavy fermion x with mass M 3> m, [i 2 /m and the same coupling to A^ 1 as if). 

The interaction preserves a global U(l) symmetry with charge +1 for (f>* and ipl and 

charge —1 for tp^- Integrating out \ we obtain 11 ! 

4,t = ^(0f 2 + O> 2 ) , h = ?£ (3.19) 
where the composite operator 

= ^1(1 + ^ = ^1 (3.20) 

(that has charge —2 under the global U(l) symmetry) will have a fractional anomalous 
dimension at low energies. 

Since this symmetry acts chirally on the fermionic field one may be concerned about 
anomalies. Indeed, the axial current j 5fl = ip^rfip has a non-zero divergence [73], and in 
the Schwinger model the symmetry transformation gives rise to an additional term in the 
Lagrangian proportional to 

C e oc d»£ = --e^A v = ~^F^ (3.21) 

r 71 Lix 



10 Co(a;) does not contribute a fractional power to the correlation functions of fermion bilincars. 

n In a 4D unparticle theory, the interaction corresponding to (|3.19p would typically be nonrenormalizable. 
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While this term is a total derivative, it cannot be eliminated in the Schwinger model due 
to instanton configurations in Euclidean space that do not decay fast enough at infinity 
and for which J d 2 x Cg does not vanish [77] . This leads to the existence of degenerate 6 
vacua [7H1 [79] that differ in the value of the symmetry-breaking condensate (O) oc e td . But 
in the Sommerfield model (unlike the Schwinger model that has gauge invariance) we can 
define a conserved axial current as 

f» = j5» + Ltf»A v (3.22) 

7T 

Then, since a (non-anomalous) symmetry cannot be broken spontaneously in two dimen- 
sions [HQ], there is no condensate. Note also that a non- vanishing condensate of an operator 
with a non-zero dimension would be in conflict with scale invariance. In the Schwinger 
model, the dimension of O at the IR fixed point is 0, but we will see that this is not the case 
in the Sommerfield model. 



4 Correlation functions of the unparticle operator 

The vertices in (13. 18[) that couple the standard model field <fi to the fermions of the Som- 
merfield model involve also the heavy field \. Since x effectively propagates over distances 
of order 1/M which are much shorter than the distances we want to consider, our effective 
interactions will involve two s coupling to the leading operator in the operator product 
expansion (OPE) [8lT 170] of a product of two tp s. In particular, we will be interested in 

Ti)*p(x 2 )ip a {xi) = c(x 2 - xx) ^■0 a (x 2 ) H (4.1) 

that defines the operator O = ip^ipi from (13.201) and its conjugate O* = iplip2- Using this 
OPE and the exact solution for the fermionic correlation functions from section 13.11 we will 
determine the correlation functions of O and O* . The coefficient function c(x2 —Xi) will not 
play any role since there is also a x propagator connected between xi and x 2 - Evaluating 
the convolution of the two at zero external momentum, that is all we need for energies much 
below M, would give a constant. 

4.1 2-point function 

Consider the fermionic 4-point function 

G« = (01^(^)^(^1)^(2/2)^(2/1)10} (4.2) 
We can express it as the free-fermion skeleton in figure HI that is 

G hl = iS o 1 ( x i - x 2) iSo 1 (Vi - 2/2) ii« 2 <W 2 - iS o 1 ( x i - 2/2) ^o 1 (Vi - X 2) ^ ai f3 2 <5/3 lQ2 (4.3) 
multiplied by the bosonic factors ( 13.101) . where 5 a /3 Sq(x) is the free-fermion propag ator@ 

= / (2^2 ^ = -^^2—- ( 4 - 4 ) 

12 We define the propagator S a (x) as 5 af3 S a {x) = -i (0\Tip a (x) ^S(0)|0>. 
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X\, 0L\ 



Xi, 0L\ 



x 2 ,a 2 yi:/3l x 2 ,a 2 yi) p 1 

Figure 4: The free fermion skeleton corresponding to the 4-point function (14. 2 j) . 

X!,a x y 2 ,(3 2 



+ 



+ + 



x 2 ,a 2 yu p 1 
Figure 5: The sign factors (77^, Kij) for the case a\ = (3 2 , a 2 = (3\ (a\ 7^ a 2 ). 



S 2 (x) 



d 2 p 



-ipx 



V 



X 



p z + xe 



2ti x 2 — ie 



where we use the lightcone coordinates 13 ! 

x ± = x° ± x 1 

For the O and O* operators, we are interested in the case 

Oil = Aj , "2 = A («1 7^ «2] 



(4.5) 



(4.6) 



(4.7) 



Only the second term of (14.31) survives, and the corresponding rj and k factors of (13 . 1 0[) are 
shown in figure [5j In the limit X\ — Y x 2 = x, yi — > y 2 = y this gives 

(4) _ C (xi - x 2 ) C (yi - y 2 ) 



On the other hand, using the OPE (14. ip . we can write (14 .2p in the same limit as 

= c(x 2 - xi) c(y 2 - yi) (0|TO(x) O*(y)|0> 
Comparing (I4.8P and (14. 9p we see that we can take c(x) = C (x)/C(x), and 

C(x) 4 



iA (x) = (0\TO(x) O*(0)|0) = C(a:) 4 iSg (ar) iS£(ar) 



(2tt) 2 (-x 2 + ie) 



(4.1 



(4.9) 



(4.10) 



13 Some properties of the lightcone coordinates are listed in appendix 151 
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At distances large compared to I/mousing (I3.17p . we see that the 2-point function is 



proportional to an unparticle propagato: 

1 

(27r) 2 (£m) 2a (-x 2 + ie 

where 



iA {x) iA u {x) = - — — — Nl+Q (4.1i; 



e 2 1 

a= . = 2T7 (4.12) 

7rrrr 1 + 7rmQ/e z 

denotes the anomalous dimension of the composite operator O = i/j^i (its total dimension 
is g?w = 1 + a). For < mo < 00, a is fractional, which leads to unparticle behavior In 
momentum spaced 

iAu(p)= , ^ (a) , - (-p*-ie) a = A{a \ r dM 2 (M 2 Y- l - (4.13) 

KF> 2(£m) 2a sin(Tra) V y J 2n^m) 2a J V ' p 2 — M 2 + ie K J 

where 

A(a) = _^n-a) = 1 (414) 

v ; 2 1 + 2 «7rr(l + a) 2 l + 2a [r(l + a)] 2 V ; 

Since 

A(a 



I-A w (p) = -^W«(P 2 ) Q (4.15) 



the unparticle phase space is 

^p) = ^ a (pr°(p°)°(p 2 ) (4-i6) 

Note that in the Schwinger model limit the phase space vanishes since A{— 1) = 0. 
Because we have the exact solution, we can write (I4.10p for arbitrary x as 

00 (_a \ n 

tA (x) = iA u (x) exp [-4maA(x)\ = iA u (x) ~ — [iA(x)} n (4.17) 

n=0 U ' 

At distances not large compared to 1/m, the higher terms in the sum in (14.171) become 
relevant. This will be important for studying the transition between the unparticle and 
particle regime. 

4.2 Higher n-point functions 

In a similar way, for a general 2n-point function of the operator O we have 

(OlTOCan) . . . 0{x n ) 0*{y x ) . . . O*(y n )\0) = \ 3 * g££ (4.18) 

1 Lk>j ° \ X j ~ X k) ° \Vj ~ Uk) 



14 Here and below, we incorporate a dimensional factor of l/(£m) 2a in the unparticle propagator so that 
it has the same engineering dimension as the O propagator. 

15 For mo = we obtain a = — 1, i.e., du = 0, which is the Schwinger model result [82] . 
16 The Fourier transform is calculated on p. 284 in Ref. |83j . 
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where is the corresponding free fermionic 4n-point functional 



r (4n) _ ("I)" IW^j ~ XkfiVj ~ Vk? 

free i^ 2 ) n 1 !,,(•'•, - ///.O 2 1 ' j 



Thus, using fOOjl . 



(0\TO( Xl ) . . . 0{x n ) 0\ Vl ) . . . O*(y n )\0) = = .^f ^ , (4.20) 

llk>j i&oixj - x k ) lAoiVj - y k ) 

The 2n-point function involves the two-point function iAo{x), and also its inverse 

tAo(x) = (4.21) 

We can represent the 2n-point functions diagrammatically with 2n vertices connected by two 
kinds of lines: solid lines representing iAo(xj — yk)', dotted lines representing iAo(xj — Xk) 
or iAo(yj — yk)- For example, the 2-point function is just a single solid line, and the 6- and 
10-point functions are shown in figure [6j 




Figure 6: The 6-point function and the 10-point function of O. 

It is important to remember that these diagrams, while useful as reminders of the struc- 
ture of (|4.20p . are not Feynman diagrams. Eq. (I4.20p is the sum of an infinite number of 
Feynman diagrams and describes the full 2n-point function^ 

Although we will not always indicate it explicitly, one should remember that the solid 
lines are actually directed, carrying the conserved chiral charge from an O* vertex to an O 

17 Note that this simple form arises from the sum over all possible contractions of the free fermion fields. 
18 This includes also diagrams that can become disconnected in certain limits, since we included such 
diagrams in (|4.19p . 
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vertex. The dotted lines are not directed. So for example we could represent the 6-point 
function as in figure [7J Curiously, the charge carried by the solid lines depends on n. In 
the 2n-point functions, each solid line carries 2/n units of charge. It is obvious that this 
conserves charge, because n lines emerge from each O* vertex and n flow into each O vertex. 





u» » » * * * n 




Ql* * ■ 



















Figure 7: The 6-point function with the direction of the charge flow indicated. 

In this paper we will do most of our analysis on the O 4-point function: 

(0\TO(x 1 )O(x 2 )O*(y 1 )O*(y 2 )\0) 

(4.22) 

_ iA (xi - 3/1 ) iA (xi - y 2 ) iA (x 2 - yi) iA (x 2 - y 2 ) 
iA (xi - x 2 ) iA (yi - y 2 ) 

While its form may look peculiar, especially because of the propagators in the denominator 
that lead to IR divergences that we will analyze in section I4.3[ it is a typical form for 
conformal field theories, in any number of spacetime dimensions. More specifically, conformal 
invariance restricts the 4-point function of any four primary scalar operators Oi with scaling 
dimensions di to have the form [49] 



(0|TOi(xi)O 2 (x 2 )O 3 (x3)O 4 (a;4)|0) = F(n,t;)n(4) dT/6 " ( * + ^ )/2 (4-23) 



where dx = di + d 2 + d% + d^ and F(u,v) is an arbitrary function of the conformal- invariant 
cross-ratios 

2 2 2 2 

u = ^f, v = ^f (4.24) 

X 13 X 24 X 14 J '23 

When all di are equal, like in our case of interest, this can be written more simply as 
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|TOi(xi)O 2 (a; 2 )O3(x3)O4(x4)|0) = J^fJ (4.25) 

\x 13 ) \x 2li ) 



In the unparticle limit of our model, i.e., with Aq — > Ay in ( I4.22p . our 4-point function 
matches this form with f(u, v) oc v d . The function f(u, v) can be different for other operators 
and in other conformal theories. However, because of the form of u and v in (14 .24 p . when 
we expand f(u,v) in powers of u and v each term will generically have powers of some x 2 j 
in the numerator and some in the denominator, similarly to what we have in our model. 



4.3 Momentum space and IR divergences 

In momentum space the 2-point function (14.101) [recall also (14.171) ] becomes 



a (-4vra)" f d 2 p u 

n=0 J V ' 



n j2 



(2n) 2 6 2 \p-p u -J2 Pj 

(4.26) 

This describes a sum of two-point diagrams in which the incoming momentum P splits 
between the unparticle propagator and n massive scalar propagators^ The discontinuity 
gives the phase space 

(4.27) 



$(P) 



A(a) 



E 



vs ' n=0 
n J2 



nl 



(27T) 



X 



n 



(2tt)^ \P-pu 




This is the form that we used in [53] to discuss the transition from the low-energy unparticle 
regime and the high-energy free-fermion regime in our toy model. We will review this in 
section [5j But to go beyond this simplest process, we need to understand the meaning of 
the higher n-point functions ( I4.20p in more detail. 

The two point function Aq(x) goes to zero as — x 2 — > oc. But that means that its inverse 
(I4.2ip goes to oo as — x 2 — > oc. At first (and perhaps second) sight, this looks like a recipe 
for infrared divergences. Indeed, the Fourier transform of Aq(x) is very singular as p 2 — > 0. 
Formally, the dotted line propagator in momentum space can be written as an expansion 
similar to ( 14T26|) : 



iA (k) 



Akx 



d 2 x 



iAo{x 



■ikx 



d 2 x ■ 



iAiAx 



■ exp [47riaA(x)] 



— i- 



.8tt 4 A(-2 - a)(£m 



,2a 



sin M 



d 2 p 



(2tt) 2 (-p 2 - ie) 2 + a 



E 

n=0 



{An a) 7 



TV. 



n 



d 2 Pi 



■iA(pi 



(27Tf5 2 (k-p-^PJ 



(4.28) 



The inverse of the unparticle propagator, which gave the 

1 



-p 2 — ie) 2+a 



(4.29) 



19 Similar behavior in the Schwinger model, where the unparticle stuff is absent, has been discussed in 
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factor in (14.281) . will lead to infrared divergences. The Fourier transform 



2-2(1+7) r (_ 7 ) 



d 2 x e*» {-x 2 + lef = _____ (4.30) 



tt r(l + 7 )7 v y (-p 2 -ie) 1+7 

has a well-behaved Kallen-Lehmann representation for — 1 < 7 < 0: 

sin(7T7) /"°° ds 1 



Jo sl+7 p 2 — s + ie 



(4.31) 



The spectral function is positive and the singularity at s = is integrable. We used this 
implicitly in ( I4.27P with 7 = — 1 — a. However, for the dotted line we needed to take 7 = 1 + a 
to obtain ( I4.28j) : this means < 7 < 1, and then the integral in (I4.3ip does not converge 
near s = 0. More practically, the factor ( 14. 29ft will lead to divergences in cross-sections. 

A better approach for doing the Fourier transform of the dotted line is to use ( I4.30p with 
7 = a and get an extra factor of x 2 by differentiating both sides as 

o-2(2+a) Y(_ n \ r 1 



tt T(l + a) J v ' dp+dp- (-p 2 - ie) 1+a 

where we use the lightcone components ( 14.61) . It looks naively as if the derivatives just 
reproduce (14.291) . Indeed, for p 2 ^ 

d 2 1 x9 1 N 

-(! + a ) 2 7—r— (4-33) 



Qp+Qp- (—p 2 — ie) 1+a (—P 2 — it 

But it would be too naive to assume ( 14.331) for all p 2 because of the singularities of these 
generalized functions at p 2 = 0. The derivative operator clearly produces a well-defined 
generalized function, but the singularities at p 2 = would require special care. 

We can think of two ways to deal with these infrared issues. In some situations, we 
can use (I4.32p directly by routing external momenta through the diagram in such a way 
that we can take the derivatives outside the loop. Then we can use (14.301) and (14.311) . The 
other alternative is to use (14.321) with an explicit infrared cut-off. This second alternative 
is instructive as well as useful, so we will outline it here. We will discuss both methods in 
more detail in section 17.11 and appendix 

Let the function f\(s) be an IR-regulated version of l/s l+a : 

for s 3> A 



/a(.s) - { sl+a (4-34) 
as s -> 



(for example, we could take f\(s) = s a /(s + A)). Then an IR-regulated version of 

7r(l + a) 2 1 
sin(7ra) (— p 2 — ie) 2+a 



(4.35) 
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IS 



d 2 



dsf x (s) 



1 



P 



+ 



s + le 
2s 



dp + dp 

^ A ^ V (p 2 — s + it) 2 ' (p 2 — s + ie) 3 



(4.36) 



s + ie 



— s + xe 

The spectral function in ( I4.36P does not have a uniform sign. In fact, it is a total derivative, 
so the integral over s vanishes. As A — > 0, the deviation of 



1 + a 



(4.37) 



from l/s 2+a is squeezed down to s = 0, but the integral over s continues to vanish. 

Because we now have a Kallen-Lehmann representation for the dotted line, we can safely 
write down the corresponding phase space from the discontinuity across the cut: 



$(fc) 



8vr 4 A(-2 - a)(£m\ 



2d 



E 

n=0 



(1 + a) 
(4vra) n 



d 2 p 



(4.38) 



TV. 



n 



d 2 pi 



■2ttS(p 2 - m 2 )9(p° 



(2-K) 2 5 2 [k-p-Y,Pi 



We will see in appendix [C] how (I4.36P and ( I4.38P work in a simple example. 



5 Disappearance process: (f> + (f> — > U 

We now review the physical process described in figure [TJi: 

+ — > Sommerfield stuff (5-1) 

Because <p 2 couples to O* at low energy (from (I3.19P ). we can obtain the total cross-section 
for this process from the discontinuity across the physical cut in the O 2-point function. 
The cross-section for a given initial state / to scatter into any possible final state F can be 
obtained from the amplitude of I going back to / using the relation 

[ d ®F \M (I -»■ F)\ 2 = - Disc iM (/ /) (5.2) 

F ^ 

where is the phase space of the final state F, and Disc refers to the discontinuity across the 
branch cut, Ai(s + ie) — Ai(s-ie), where s = E 2 m . This is also known as the optical theorem, 
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and the right-hand side of (15. 2\i is often written as 2 Im M. (I — > I). In our particular process, 
for (j) momenta Pi and P 2 , this gives the cross-section 



lmM(P 1 ,P 2 ^Pi,P 2 ) h 2 
a = = Im A (P) 

s s 



(5.3) 



where P = P\ + P 2 and s = P 2 . In the unparticle limit (-^/i <C m), using (I4.15p . or di- 
rectly the phase space (I4.16p . we find the fractional power behavior expected with unparticle 
production: 

= A{a) h 2 1 
° 2 ^m) 2a s 1 - a 1 ' ' 

On the other hand, in the free-particle limit that appears at high energies \fs 3> m, we have 
C(x) -+ 1 in 0| and then 

h 2 1 

(5.5) 



(7 



4 S 



which is the cross-section for + — > ip 2 + ^1 without the Sommerfield interaction. 
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Figure 8: The cross-section for the disappearance process, based on the expansion (j4.27p . The long- 
dashed lines are the standard model particles 4>, the solid line represents the unparticle propagator while 
the short-dashed lines represent massive boson propagators. The thick line is a cut for computing the 
contribution to the cross-section. 



The transition between the two limits can be studied using the expansion in (I4.27P which 
shows that the general result can be described as the production of unparticle stuff along with 
an arbitrary number n of massive bosons (to the extent that this is allowed energetically), 
see figure [HJ Along with an additional massive boson, each subsequent term is proportional 
to an additional power of a. One can easily obtain explicit results in the case of small a, 
when only the first few terms in the expansion contribute. For example, the n = 1 term 
contributes 

$« = -ae(^S-m)\Yi^ + 0(a 2 ) (5.6) 

m 

which gives the total phase space at the leading non-trivial order in a as 



1 



— a 



In 



2 



j + 6{\fs-m) In 



m 



+ 0{a 2 ) 



(5.7) 
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For energies y/s > m, this expression reduces to 

$ = I + 0(a 2 ) (5.8) 

that is the free-fermion result ( 15. 5p . Thus, for \a\ <C 1 there is a discontinuity in d^/dy/s at 
y/s = m, where a transition occurs from pure unparticle behavior below energy m to pure 
free-fermion behavior above m (see figure [9]) @ Interestingly, the free-fermion behavior is 
obtained here as a sum of the n = and n = 1 terms of figure [HJ 



0.5- 




Figure 9: Phase space for the disappearance process in figure [Tk as a function of the energy yfs (in 
units of m) for a = —0.1. 

For larger values of \a\, higher powers of a will need to be taken into account in the 
contributions with extra bosons in order to reach the free-fermion regime at sufficiently high 
energies. Since each additional massive boson that we include gives a contribution with one 
extra power of a, we will need to include N massive bosons to cancel all a-dependent terms 
up to 0(a N ). Then the free-fermion behavior will only appear at y/s > Nm (while in the 
range m < ^/s < Nm mixed behavior will be observed, with discontinuities in d<&/d\fs for 
each multiple of m). In the limit a — > — 1, the required value of iV becomes infinitely large: 
the condition \a | 1 implies iV > _ ln |_ a - ) ■ Note also that for a = — 1 the unparticle 
contribution disappears since A{— 1) = 0. This limit is the Schwinger model; it has been 
studied in [H|82]. 

6 Missing charge process: + 0— > (j) + (j) + U 
6.1 Inclusive treatment: general 

In this section we show that it is possible to go beyond the simple calculation of section [5] 
to study processes in which standard model particles are radiated from the unparticle stuff. 
We begin with a process in which infrared issues do not intrude at low energies. In section [7] 

20 The linear approximation (|5.7[) is not valid for ^/s <C m due to large In -^/s, but we have the exact 
expression (|4.16[) . 
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we will consider another process in which the IR properties of the dotted line are more 
immediately important. 




Figure 10: The 4-point function (I4.22|) , in position space on the left and momentum space on the 
right. Momenta k and k' enter the diagram and k and k leave it. The thick line is a cut for computing 
the cross-section of the process (16.11) . 

Consider the missing charge process of figure [lb, 

(f>{pi) + 4>{P2) ~> 0(9i) + 0(^2) + Sommerfield stuff (6.1) 

Because of (I3.19p . and the fact that four units of U(l) charge must be carried by the Som- 
merfield stuff in this process, the Sommerfield dynamics that contributes is associated with 
the discontinuity across the cut in the O 4-point function (I4.22p that is represented by fig- 
ure [TOj We obtain one contribution to the total cross-section of (16. ip by annihilating two s 
at yi (and injecting k — pi + p 2 ) and creating two s at y 2 (and injecting k! = —q\ — q 2 ), 
and analogously for x\ and x 2 . More generally, k can be taken to be the sum of any two 
out of pi, P2, —qi, —q 2 , and k' the sum of the remaining two. Independently, k and k can 
also be each assigned two of these four momenta. All the possible combinations are summed 
over in the cross-section. (In other words, we are taking into account the interference of all 
the possible ways of attaching the standard model particles.) Let's define the cross-section 
factor a as 

a^^TlMI 2 ^ (6.2) 

i 

where A4i are the amplitudes of the various processes that contribute to our inclusive cal- 
culation and $j are the phase spaces of the Sommerfield stuff that is produced in these 
processes. In order to obtain the actual cross-section one would still need to multiply a by 
the couplings to the standard model h A , divide it by the usual factors involving the momenta 
of the incoming particles, multiply by the phase space factors of the outgoing particles and 
integrate over the phase space. Based on (15.21) . we have 

d 2 h dH 2 x d 2 £ 3 

0~ = y I <P(^ -l X - l 2 - «3j n^ljTTT^ $ ^2 TTT^T $ (4 



% i3 



(2tt) 2 

x iAo(h £ 2 ) iAofc -£1- h) (6.3) 
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where Ao{k) is the dotted line propagator in momentum space given by (I4.28P (with appro- 
priate IR caveats), $(p) is the phase space of the Sommerfield stuff (along with its massive 
bosons) as given in (14.271) . the momenta fc, and kj take the values 

ki, kj e {pi + p 2 , -qi - q 2 , pi - qi , pi - q% , vi - q\ , vi - qi\ (6.4) 

and we define 

P = p 1 +p 2 , q = qi + q 2 , Q = P-q (6.5) 
In the unparticle limit (i.e., without the massive bosons), ( 16.31) reduces to 



1 

a 



2 3 (2tt) 



A{a) 
(£m) 2a 



POO 

V / d£+ d£+ d£+ d£± d£^ d£ 3 



x9((Q-£ 1 -£ 2 - £ 3 ) + ) 9 ((Q -£ 1 -£ 2 - e 3 )~) 

x WPi{Q - h - 4 - hYK^W -h-h- hY\ a 

x Aoih - h - i 2 ) Aoikj - li - h) (6.6) 



6.2 Results for small a 

It is easy to compute (16.61) for \a\ <C 1, since then 



^o(k) = - 2 + 0(a 2 ) (6.7) 



and 



a 



i r°° 

^J2j o d£+ditd£id£^dqdq 



2 7 (2tt) 



x9({Q-£ x -£ 2 - ((Q 

x A (k { - - £ 2 ) Aoikj - £, - £ 3 ) 



2 



where^l 

/>oo 

I) = y d£ 1 d£ 2 d£ 3 ^ f _ 



oo 



o (A; - — 4) 2 (k — £i~ 4) 

(Q - 2jfe) In (1 - Q/k) - (Q - 2k) In (1 - Q/k) 
(k-k){Q -k-k) 



9{Q) (6.9) 



21 The following integral would be divergent for < ki < Q or < kj < Q, but the kinematics ensures 
that ki and kj never fall in these ranges. 
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To analyze (I6.8P in more detail, we can consider, for simplicity, only the termini 

k h kj E {pi + p 2 , -qi - q2} (6.10) 

while omitting all the other interference terms in ( 16. 4p . Notice that for terms with k = k or 
k = Q — k, (I6.9p reduces to 



I{Q,k,k) 



Q 21n(l-Q/Jfe) 



9{Q) (6.11) 



--7^1 (6-14) 



_k(k-Q) 2k -Q 
which gives, with the definitions (16. 5p . 

9«/^m fP + ~q + 21n(<7+/P+)\ (P- -q- 2In(qr/P-)\ 

For q ^ P, that is when most of the energy goes into the unparticle stuff (which is possible 
if P ^> m^), (7 is dominated by the term 

aoc^ (6.13) 

This small-g enhancement occurs because the momentum flowing through the dotted line 
(whose propagator is iA {k) ~ \/{k 2 ) 2 ) is allowed to be as small as q. On the other hand, in 
the limit when only a small fraction of the momentum goes into the unparticle stuff, 
the cross-section behaves as 

- „ 

(P 2 

This result can be understood by observing that in this regime the points x\ and X2 in 
figure [10] are typically very separated from y 1 and y 2 (compared to the separation between 
x\ and X2 or y\ and y?). Then the four solid lines are essentially connected between the same 
two points and given by [zAjy(X)] 4 oc (— X 2 )~ 4 ( 1+a ) (where X = x — y). In momentum space 
this becomes oc (— Q 2 ) 3+4a . Each of the dotted lines has iK u {P) oc (— P 2 ) _2_a , so we expect 

° « (p2)4+2a ( 6 - 15 ) 

This indeed agrees with (I6.14p (in the limit |a| <§C 1 that we assumed in the derivation of 
(16.141) ) . We will discuss this point of view in more generality in subsection 16.31 

The missing charge process in the unparticle limit does not require us to use the IR- 
regulated form of the dotted line because the kinematics keeps the momentum carried by 
the dotted lines away from the light cone. At higher energies, above threshold for the 
production of the massive bosons, we do have to worry about these IR issues because while 
the total momentum carried by the dotted line cannot be lightlike, the momentum carried 
by its unparticle part can be. We will not discuss this here, because we will see a related 
issue already at low energy in the process analyzed in section [71 



22 In particular, these would be the only terms present if instead of (|3.19p we considered a model with 
two flavors of the standard model scalars, (j>A and 4>b, with C- m t = \ \P (4>*a 2 + 4>b 2 ) + ^* {'Pa + 0fl)] an( l 
asked about the process 4>a{pi) + <Pa(p2) — > ( / ) b( ( Zi) + 4>b{ < 12) + unparticle stuff. 
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6.3 Inclusive treatment: series expansion 



k ^ k 

Q 



Q 



k 



■Py- 



k' 



Figure 11: A single term in f)6. 17[) . The p x and p y lines are described by (j6.18p and the Q line by 
(I6.2ip . The arrows indicate the conventional direction of positive momentum. The thick line is a cut 
for computing the contribution to the cross-section. 

We can also obtain analytic results, in the unparticle limit, for an arbitrary a, by expand- 
ing around the limit where the fraction of the momentum that goes into the unparticle stuff 
is small: Q <C P. In this limit the contribution should come mainly from configurations of 
figure [10] in which ( x ,(y X, where 

( x = x 1 -x 2 , C y = y 1 -y 2 , X = x 2 -y 2 (6.16) 

Expanding the 4-point function (14. 22ft in ( X /X and ( y /X we have 

(0\TO(x 1 )O(x 2 )O*(y 1 )O*(y 2 )\0) = iK u {Q i~^u{Q [i^u{X)f 



x 



2(1 



+(l + a 



,C Cy +4(1 + g)2 (Co ; Cy? 



x- 



X 2 



(3 + 2a) (C 2 + C,- 2 ) - (5 + 4a) C x C y 



(x- 



, 2a\, 2(2 + a)(c 2 + C 2 )-(5 + 4a)CC, 
-(l + a)(l + -J (C-Cy) — 



-2(1 + a) 2 (C-C) 



(X-) 3 

(3 + 2a) (C 2 + C y 2 ) - (5 + 4a) ( x ( y 



{x-yx 

+ {— -H- + for all asymmetric terms} + . . . 



(6.17) 



Each term in (16. 17p describes a diagram as in figure [UJ in which momentum Q flows 
through X, momentum p x = k through ( x , and momentum p y = —k through ( y . To 
take the ^-dependent factors into momentum space, note that the Fourier transform of 
*A w (C)(C) m (CT ^ 

m— 1 n— 1 



'2 d 



i dp A 



2 d 



i dp~ 



iA u (p) = (2t) m+n Y[(2 + a + k) ]J(2 + a + k) 



k=0 



k=0 



(p+) m ( P -y 
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(2i) 



m+n 



T{2 + a + m)T{2 + a + n) ih u {p) 



where, as in (14.2811 . 



iAw(p) 



[r(2 + a)] 2 
2 3+2a (27r) 3 r(2 + a)(^m) 2a 



(p+) m (p- 



(6.18) 



(6.19) 



r(-l - a) (-p 2 - it) 2+a 

To obtain the contribution of a term from (16. 17ft to the cross-section, we take the product 
of two factors of the form (16.18P coming from the ( x and ( y lines and multiply it by the 
discontinuity across the cut through X. To find the latter, we compute the Fourier transform 



(X~) M (X+) N 
(suppose, for definiteness, that M > N) as 

\M—N\a A CVM4 /« Q \ M—N 



(6.20) 



/ 



(X+)^" w [zA w (X)]^ 



(X 2 - * e ; 



2 9 



z <9Q- 



d 2 X 



e'Q x [zA^(X)] z 
(X 2 - ie) M 



(-1) A H 



r(-3-4a-M) /2 9 



2 7+8a+2M( 27r )7(f m )8a f(4 + 4a + M) 0Q 



M-N 



(-Q 2 - ie) 



3+4a+M 



;l+M+iV 



T(-3-4a-X) 



(Q + ) M - N (-Q 2 -ie) 



2 7 + 8a+M + N( 27r y(£m) Sa T(4 + 4a + M) 
which gives the phase space factor 



3+4a+V 



(6.21) 



;M+Nt 



-1) N sin(47ra) T(-3 - 4a - N) 
2 6 +8a+M +iV(27r) 7 (£m) 8a T(4 + 4a + M) 



(Q + ) M (Q-) 1 " (Q 2 Y +/ta e(Q°)9(Q 



, N 



2 \ 3+4a 



:M+N 



(Q + ) M (Q-f (Q 2 f +/la 0(Q°)9(Q' 



2 7+8a+Af+iV(2 vr )6^ m )8a f ( 4 + 4fl + M ) p( 4 + 4a + jy) 

Summing the contributions of all the terms in (I6.17P and adding an overall minus sign from 
(15.2p . we obtain 

A u (p x )A u (p y ) (Q 2 ) 3+ " a 0(Q°)0(Q 2 ) 



a 



x 



2 7+8a ( 27 r) 6 (£m) 8a [r(4 + 4a)] 2 
1 1 



M Py 

(2 + a) 2 Q+ 2 (2 + a) 2 



(2 + a)(3 + a)(3 + 2a) 



4(5 + 4a) 



+ 2 



+ 



P. 



PIP 



+ 



fx try 



1 



(2 + a) 2 (3 + a) 
24(5 + 4a) 

(2 + a) 2 



" 1 vi J 



2(4 + a 



1 1 



Pf 



Pf 



1 + 4a 

P^ 



Q +2 Q- I — - — 

Px 



1 



1 



(3 + a)(3 + 2a) 



5 + 4a 

{+ <H- — for all asymmetric terms} + . . . 



Pf 



Pf 



2 + a 

PtPy 



(6.23) 
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Summing over the 4 possibilities in (I6.10p we get 

lK u {P)?{Q 2 ? +Aa e{Q Q )em 



a 



2 5+8a (2vr)6(£ m )8a [r(4 + 4a)] 2 

X 



Q + , Q~\ , , n , ^Q 2 



1 + (2 + ( + ^ir ) + (2 + 



+ 



(2 + a)(28 + 31a + 8a 2 ) 
4(5 + 4a) 

(2 + a) 2 (3 + a)(17 + 8a) 



cr 



-\ 2 



+ ( 2 + a ) ( -p^ + pip- 



er 

P^ 



12(5 + 4a) 

Note that at the leading order in a this becomes 



Q_ 

P- 



+ 



(6.24) 
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a 2 (Q 



2\3 



18 (ps 



+ 
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23 / q 2 q 



51 
p- 



O 2 14 

+ 4— H 

P 2 5 



0^ 
p+ 
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p- 



Q 2 Q- 

5 lP 2 P+ ' P 2 P- 
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P^ 



+ 
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p- 



+ 



(6.25) 
0(Q°)9(Q 2 ) 



which is precisely ( I6.12p expanded in powers of Q/P. 



6.4 Interpretation in terms of exclusive processes 

What physical states of the conformal sector are produced in the missing charge process 
whose inclusive cross-section we have computed here? Can the inclusive sum in (16. 2p be 
decomposed into distinct well-defined contributions? 

To understand the answer to these questions, let's look at our calculation from a different 
point of view. Consider the series expansion that we did in subsection 16.31 Since we took 
the limit ( x ,(y ^ X, we can replace the products 0(xi)0(x2) and O* (yi)0* (2/2) in fl6. 17j) 
by their OPEs. In particular, using the shorthand notation 



(a\b) = (0\Ta{X) 6(0)|0) 

we can write (I6.17P as 

(0\TO(x 1 )O(x 2 )O*(y 1 )O*(y 2 )\0) 
ih u {( x )ih u {Q 
(0 2 \0* 2 ) + ~ [Q(d„0 2 \0* 2 ) + Q{0 2 \d^O* 2 }] 



(6.26) 



3 + 2a 
4(5 + 4a) 



Q 2 (d 2 _ &\o* 2 ) + c\&\d l _ O*' 2 ) 



- 2 //n2i q2 sn*2\ 



~ A QQ{d-0 2 \d„o* 2 ) 



+t [Cc;(d + o 2 \d^o* 2 } + cc(d + d„o 2 \o* 2 } + c+c y (o 2 \d + d_o* 2 }} 
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2 + a 
+ 12(5 + 4a) 

3 + 2a 
+ 8(5 + 4a) 



C~\dlo 2 \o* 2 ) + c~"(c 2 |^c* 2 ) 



- 3 //n2i n3 /n*2\ 



+Cc 2 (9 + 9!o 2 |o* 2 ) + c 2 C(0 2 |5 + 9!o* 2 ) 



+g [CCC;(9 + ^o 2 |9_o* 2 ) + cCC(^o 2 |9 + 9_o* 2 )] 

+ {— -H- + for all asymmetric terms} + . . . 
where we defined O 2 to be the leading operator in the OPE for 0(xi)0(x 2 ) with 



(0 2 \0* 2 ) = [iAu(X)\ 4 



(6.27) 



(6.28) 



and the leading terms in the OPE are 
TO«)O(0) = iA(C) 



+ 2 ^2 



C 3 d 3 _ + C +3 ^ 



+1^0 9 + 2 + a 
+ 4 C ^- + 12(5 + 4a) 

+ J^tC 2 fc d • c/,7.)^7 • ... 



(6.29) 



8(5 + 4a) 



C 2 (0) + ... 



The coefficients of the derivative terms in (I6.29|) must in fact have this precise form in 
order for the 3-point function {0\TO(x)O(y)O* 2 (z)\0) to be consistent with the conformal 
symmetry. In general [HS] (see also [56] EH EZ]), the contribution of all the derivatives of a 
scalar operator of dimension A (in our case, the operator O 2 has dimension A = 4(1 + a)) 
to the OPE of two scalar operators of dimension d (in our case d = 1 + a), in D spacetime 
dimensions, should be proportional to 



(C 2 ) 



is/5 f dt - *)] A/2 ~ X ofi (a + 1 - f ; (i - C 2 □) 

r(A+ 1 -D/2) (-C 2 n) m 



((2)^-A/2 Z_, 4 m m! F ( A + X _ D /2 +m ) y Q 

1 ? i°+^ + ^±4,(c^) 2 - 



dt[t{l-t)^ /2+m - l e^ a 



(C 



2 xd-A/2 
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4 + A 



48 (1 + A) 



(1 + A) 
(C • df 



16(1 + A)(l + A- D/2) 
A 



C 2 D 



32 (1 + A) (1 + A - D/2) 



C(C-d)D + 



(6.30) 



Furthermore, there exists an explicit expression for the contribution that each rank-£ 
tensor Ok in (12. 3 p (together with all its derivatives) can make to the 4-point function of 4 
scalars 1601 ■ The contribution to 



(0\TO 1 (x 1 )O 2 (x 2 )O 1 (x 3 )O 2 (x i )\0) 



(6.31) 
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where 0\ and O2 have the same dimension d is proportional tcF*l 

+ {^^r/ + } (6.32) 
where 2 -Pi is the ordinary hypergeometric function, 

X 13 J '24 

and u = rj 2 as in (I4.24p . In our case 

a?i2 = Cx , ^34 = C y , x 2 4 = X , xi 3 = X + - Cy (6-34) 

Expanding H 6 . 3 2 1) in small ( x and C^, with i = and A = 4(1 + a), we obtain precisely the 
expression (I6.17P (up to an overall prefactor that is not fixed by the conformal symmetry). 
No other primary operators besides O 2 are required for reproducing the terms shown in 
(16.17p . and this confirms the interpretation (16.271) . It is then obvious from ( 16. 2 7ft that taking 
the cut in the 4-point function as we did describes the production of O 2 stuff: 

+ ^0 + + {O 2 stuff} (6.35) 

Additional operators do appear in higher-order terms that are not written in (I6.17p . In 
particular, at order ((/X) 4 there appears the extra contribution 




iA u (QiA u (( y ) VA U {X)Y A__±_ ) (( „ )(i) 

that must be accounted for by the two-point functions of some new operator(s) of dimension 
A = 4(1 + a) + 2 and rank t = 20 

In conformal theories that are not solvable, which is typically the case in four dimensions, 
there is no easy way to determine which operators contribute to the OPE. However, the 
conserved currents of the CFT are the usual suspects to appear in OPEs and they always 
have their canonical dimensions. In our model, they will indeed appear in the OPE for 
G*(C) O(0) that will be relevant for the missing energy process in section [71 The conserved 
current in our model could not appear in the OPE for 0{Q O(0) because O is charged 
under the axial U(l) symmetry while j M is not. Yet further information can be obtained 
based on the conformal symmetry alone. For example, an upper bound on the dimension 
of the leading scalar operator in the OPE of two identical real scalars was derived in [88J. 
In our case, considering the operator O + O* (we are taking this combination in order to 



23 We present the expression that is relevant to 2 spacetime dimensions, but an analogous expression for 4 
dimensions is given in [60] as well (see also the explanations in [88]). 

24 In fact, it can be shown by more advanced methods that all operators of vector charge and axial charge 
4 in this theory have dimensions of the form A = 4(l + a)+n where n is an integer [55] , 
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have a real operator) this bound implies that the OPEs for 0(C) O(0) or 0*(C) O(0) must 
contain a scalar operator with dimension below 0.53 + 4(1 + a). This condition is satisfied by 
the operator O 2 (whose dimension is 4(1 + a)) that appears in 0(C) O(0). For a > — |, the 
condition is satisfied even earlier by the dimension-2 operator j 2 that appears in 0*(C) O(0) 
(we will analyze this OPE in section ITU]) . 



6.5 Exclusive treatment: amputated 3-point functions 



We can also compute the exclusive cross-section for (16.351) directly using the amputated 
3-point function as discussed in section [2j For D = 2 and d = 1 + a, (12 . 1 9[) reduces to 



|M| 2 $ 
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Wick rotating fc 1 , the integral over becomes 

/c dk 2ni 
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and then 
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(6.39) 



(6.40) 
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The integral in (I6.40p is hard to solve analytically. However, expanding the integrand in pow- 
ers of Q/P it can be computed [90] (assuming for the purpose of calculation the unphysical 
range —2a < A < 0) to give 

r(l - A) [r(A/2)] 2 T(A/2 + a) 1 
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so 
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(6.42) 



For the process (16.351) we have 

A = 4(1 + a) , C 2 



C, 



1 



(27r) 8 (£m) 8a ' J (27r) 6 (£m) 6a 

where A and C 2 follow from (16. 28 p . and C3 from (I6.29p and conformal invariance. Then 
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(6.44) 



After multiplying this by a factor of 4 to account for the various ways of attaching the 
standard model particles (I6.10p . this agrees with (I6.24p . 
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7 Missing energy process: + 0— > (p + (p + U 
7.1 Inclusive treatment: general 



0(x 2 ) 0*(y 2 ) 



0*( Vl ) 0{x x ) 




P-q- 



P-q 



1—12 — t-3 



'-I — *2 — ^3 



Figure 12: The 4-point function (I4.22D , oriented in a way relevant to the process (17. 1ft . Momenta P 
and — g enter the diagram at the bottom and leave it at the top. The thick line is a cut for computing 
the cross-section. 



Consider the missing energy process of figure [£b, 

+ 0G°2) — > <P(qi) + 0(92) + Sommerfleld stuff (7.1) 

Because of (13.191) . the Sommerfleld dynamics that contributes is associated with the dis- 
continuity across the cut in the O 4-point function (14.221) that we now orient as shown in 
figure [T2l where we again use the definitions (I6.5p . The resulting cross-section factor (defined 
as in (16.21) ) can be computed by (15.21) as 




xiA (P -h-t 2 )i^o{P -Z1-I3) (7.2) 

where Q(k) is the Sommerfleld stuff phase space (14.271) and $(fc) is the "phase space" corre- 
sponding to the dotted line (14. 38[) . 

In this process it is immediately clear that we have to worry about the IR dynamics 
of the dotted lines because the momentum carried by them can be light-like. We will not 
perform the full computation in practice, but will explain how it can be done in principle 
(numerically) using the procedures discussed in section 14.31 and exemplified in appendix ICl 
However, in subsection 17.21 we will be able to obtain detailed analytic results by expanding 
around the limit of small unparticle momentum. 

Let us argue, in two different ways, that (17. 21) gives a finite result in the limit that the 
IR cut-off A in (14.381) goes to zero. For simplicity, consider the low-energy unparticle limit 
where we can ignore the massive bosons and use the unparticle propagator (14. 13[) for iAq, 
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and the phase spaces 



^i) = A^L e {t°)e{e 2 ) (e 2 ) a (7.3) 



(£m) 

16vr 4 A(-2 - a)(£m^ 2 ' 
[1 + a 

Each of the integration measures in (17.21) can be written as 



m = \ , ; w w to 2 ) + ? rm (7.4) 



= ^ (7.5) 

where s = C 2 . Now the strategy is to do the integration over the spatial momenta £j with 
the Sj held fixed. For fixed Sj, this is just like computing the cross-section for producing 
four particles with masses-squared Sj with the amplitude-squared \Ai\ 2 oc h 4 Ao(P — t\ — 
i2)A [P — £1 — £ 3 ). Because this \Ai\ 2 is a smooth function of the invariants, we expect no 
rapid dependence of the cross-section as a function of the Sj. The point is that while the 
phase space in (17.4ft is singular at t 2 = 0, the integral over s is finite and well-defined. This 
remains true when the phase space is multiplied by the smooth matrix element. 

There is also another way to see that our strategy should work and that we can remove 
the IR cut-off at the end, because there is another way to do the calculation. In figure [131 
we show a slightly modified diagram in which an additional momentum k is flowing in such 
a way that the momenta P and q flow only through the two dotted lines. We know from 
the discussion in section S3! that the appropriate definition of the dotted "propagator" is as 
a derivative. For this momentum routing, we can do the differentiations with respect to P 
and q and take the derivatives outside the integrals over the loop momenta, and then none 
of the lines requires an IR cut-off, so the result is clearly independent of A. After doing the 
loop integrations and doing the differentiations, we can set 

k = P + q (7.6) 

to obtain the result. 

One may wonder what the result of this calculation will look like. When we took the 
limit Q <C P in the missing charge process (section [6]) we saw that the phase space reduced 
to four copies of the unparticle phase space (14. 161) corresponding to the operator O. Naively, 
this could be attributed to the fact that we took a cut across four solid lines. However, this 
cannot be interpreted as the production of four units of O unparticle stuff because the O 
stuff in the disappearance process of section [5]has charge 2, while the putative "four objects" 
in the missing charge process carry a total charge of 4, not 8. Instead, it must be interpreted 
as the production of O 2 stuff (and stuff corresponding to higher-dimension primary operators 
in the OPE), that has no relation to O stuff. Similarly, the missing energy process that we 
discuss here cannot be described as the production of two units of stuff corresponding to 
the solid lines iAq(x) and two units corresponding to the dotted lines iAq(x) = l/[iAo(x)]. 
The latter would be unsatisfactory also because the "propagator" iAo(x) would describe an 
operator with dimension d = —1 — a < which is unacceptable from the physical point of 
view. In the next subsections we will see what actually happens. 
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Figure 13: Same as figure [12] but with a different routing of momentum and an additional external 
momentum k flowing through the diagram. 

7.2 Inclusive treatment: series expansion 



p - q 
\p-Q 
Q 



Q 



p 



\ 



p 



-q 



Figure 14: A single term in (17.8ft . The P and P — Q lines are described by (j7.9j) while the cut through 
the Q line gives the phase space factor (17.18ft . 

Similarly to what we did in subsection 16.31 we can easily obtain analytic results by 
expanding around the limit where the fraction of the momentum that goes into the unparticle 
stuff is small: Q ^ P. In this limit the contribution should come mainly from configurations 
of figure fT2l in which points Xi, yj with i = j are much closer to each other than points with 
% ^ j. Denoting 

d = yi-Xi, X = x 2 - xi (7.7) 
we expand the 4-point function ( 14.221) in small Q/X: 



(0\TO(x 1 )O(x 2 )O*(y 1 )O*(y 2 )\0) 
= ?A W (G)*A W (C 2 



1 + (1 + a)4%r + (1 + ^ CfCa (C " 



+(l + a) 



{x-f " ' {x-f 

crc 2 "(2cr 2 -(2-a)crc 2 " + 2C2 

2(X-) 4 
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,2 CiCf , /-i , \2C1CKC2 Ci 



+d+«)^ + (! + «) ( A - m+)2 

. . x cr 4 c 2 " - crc 2 ~ 4 -(!-«) cr 2 c 2 ~ 2 (cr - c 2 ~) 

+(! + «) 

+ { )■ + for all asymmetric terms} + . . . (7-8) 

Written in this way, each term describes a diagram as in figure [TH in which momentum Q 
flows through X, momentum — P through (1, and momentum — (P — Q) through ( 2 - For the 
(■-dependent factors, note that the Fourier transform of iA u (()((~) m (( + ) n is 

n<~i\wi/n<~i\" /<->\ m+n m—1 n—1 .» / \ 

2 d \ [2 d \ (2\ -i-r, mTT/ m «A w (p) 



i dp + J \i dp 



(n\ m+n _ _ -a / \ 

f) n(-*)ii(«-*)^^i (") 
/c — k — 

2\ m+n [r(l + a)] 2 zA u (p) 



ij T(l + a - m)T(l + a - n) (p+) m (p-) n 

To obtain the contribution of a term in (17.81) to the cross-section we take the product of two 
factors of the form (17.91) coming from the (1 and (2 lines and multiply it by the discontinuity 
across the cut through X. Apart from the first term that describes a disconnected diagram] 25 ! 
all the terms in the expansion have the form 

{x-) M {x+) N (T ' 10) 

where M and N are integers. They describe the production of some massless stuff@ For 
finding the discontinuities across the cut, consider, for example, the 1/(X~) 2 term. We 
compute its Fourier transform as 

ft e ^ (* + > 2 (I 8 V [,fX e — .' W+)2 (711) 



where we used 



\i m [d 2 X- f N9 x = -z-lim S 7 ^ (-Q 2 -^) 1 ^ 

= -^(-g 2 -.e)ln(^^) (7.12) 



25 As mentioned before, the expression we use, (|4.19l) . does include diagrams that are disconnected in limits 
where the C(x) factors in (|4.18l) do not connect them. Here the solid and dotted lines that connect the top 
and the bottom parts of the diagram cancel each other in position space in the limit y\ — > xi, 2/2 — > xi- We 
will mention the disconnected term again later in this section. 

26 Kinematically this looks like M right- moving and N left- moving free massless fermions, see (I4.4p - (I4.5|) . 
But since the M "particles" are massless and collinear (and similarly the N "particles" ) , their exact number is 
uncertain. Furthermore, some of the 1/ factors come from the momentum dependence of the interactions 
rather than describe an additional particle (we have already seen this happening in the missing charge 
process in section O . Even though terms in which both M and iV are non-zero will describe states with 
mass Q 2 > 0, there will be no dependence on a fractional power like we had for the unparticle stuff discussed 
in the previous sections. 



38 



where c is a ((^-dependent) number. This gives the phase space factor 

-2n 2 9(Q + )(Q + ) 2 5(Q 2 ) (7.13) 
In the same way, we can analyze terms (I7.10P with arbitrary M and N. Similarly to ( I7.12p . 

. .w-i n / —Q 2 — i< 



JQX 



in 



(X 2 -ie) M ~ (2 M -!(M- 1)!) : 



(-Q 2 -ie) In 



cm 



(7.14) 



where cm is a regulator-dependent number, so the Fourier transform of terms with N = 
(and analogously terms with M = 0) is 



2 9 



M 



JQX 



d 2 X 



l M + 1 IT (Q + ) 



(X 2 - ze) M 1 M -\M - 1)! Q 2 + ie 



(7.15) 



and the resulting phase space factor is 



Z M 7T 2 



(Q + ) M 8(Q + ) S(Q 2 



2 M ~ 3 (M- 1)! 

For terms with M > AT > (and similarly iV > M > 0), the Fourier transform is 



(7.16) 



2 9 



i dQ- 



M-N 



d 2 X 



% M + N+l n 



(X 2 - ie) M 

M-l /•„_\JV-1 



2 M + N ~ 2 {M- 1)\{N- 1)! 



(Q + y (qt ^ 



-Q 2 -*A (Q + ) iU (g 

T UM,N 



cm 



(7.17) 



Q 2 -te 



where the term with the unimportant prefactor d,M,N will n °t contribute to the imaginary 
part. The resulting phase space factor i^| 



i M+N n 2 



-l)\(N-l)\ {Q+)M " W")*" 1 ^^ 



2M+n-3(M - 1)\(N - 1) 



(7.1* 



Summing the contributions of all the terms in (17. 8ft and adding an overall minus sign from 
(15. 2p . we obtain 



a = 87rV(l + a) [A U (P)Y 

2 

+ (l-a) 



p+ 



+(l-a) 1 



+ {+^-} 



cr 



+ (1 - a) 1 



1 _ a(3 + a) \ /Q 



o(7 + a) \ /Q 



12 



+ • 



9(Q + )5(Q 2 ) 



(7.19) 



+ 07T a (I + a) 



p2 p2 



l + d-a)(^ + |-)+... 



e(g°) #(Q 2 
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Figure 15: The amplitude (a) and the cross-section (b) of the process (j7.20p . The solid lines are the 
unparticle propagators. 



Nicely, the IR divergences do not appear in this approach. 

We would like to note that while the cross-section (I7.19P vanishes in the limit of either 
a — > or Q — > 0, there is also another process, without missing energy, whose amplitude 
contains just the unparticle 2-point function as shown in figure [T3k : 

+ 0^0 + (7.20) 

with 

a = 2(2n) 2 [A u (P)] 2 6(Q°)5 2 (Q) (7.21) 

Figure [T5b shows the diagram that is related to the cross-section of this process by the 
optical theorem. In fact, the 4-point function knows about this diagram: it is described by 
the disconnected term in (17. 8p . The contribution (I7.2ip is hiding behind the IR divergences 
in section I7.1[ similarly to the example analyzed in appendix O 




/ \ / \ 

/ \ / \ 

/ \ / \ 



Figure 16: (a) The amplitude for the missing energy processes in the free-fermion limit, eq. (I7.22p . 
The solid lines are ifi. (b) Similarly, the free-fermion limit of the process of figure [l5b. 

In the case a = 0, Sommerfield model reduces to a theory of free fermions and the 
interaction (I3.19P describes the missing energy processes 

+ 0^0 + + Ipl + ^x, + 0^0 + + ^2+^2 (7-22) 

27 Alternatively, notice that (|7.10p can be represented as the a — > — 1 limit of (|6.20[) with £to = 2tt, and 
then reduces to (jTTg)) . 
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as shown in figure [T6k (while the process (17.201) reduces to figure [T6b). After integrating 
over the phase space of the fermions, for the first process we get 



\M\ t = -^—^ 

1 2(P+) 2 1-Q+/P 



9(Q + )5(Q' 



(7.23) 



and the second process of (I7.22p is described by the same expression with {H — > —}. This 
is consistent with (I7.19P which for a = reduces to 



_ (Q 

a - 



2(P+Y 



p+ 



p+ 



+ . . . 



(Q + ) $ {Q 2 ) + {+ -> -} (7.24) 



On the other hand, the missing charge process from section [6] does not have a free-theory 
counterpart, which is consistent with the fact that the cross-section (16.81) vanishes for a = 0. 



7.3 Interpretation in terms of exclusive processes 



We can interpret the expansion (17. 8p that we did in the Q <C X limit by replacing the 
products 0*(yi) 0(x\) and 0*{y 2 ) 0(x 2 ) in the 4-point function by their OPEs 



TO*( Vi ) 0{ Xi ) = iA u (Ci) + J> fc (Ci) Ofa) 



(7.25) 



It can be shown that in the low-energy effective theory describing the unparticle limit, the 
leading operators Ok{x) that appear in this OPE are the components of the current j 7 *, 
products of several such operators (which include the stress-energy tensor T^), and their 
derivatives (for more details, see appendix lAl where this current is denoted j^). For example, 
the 1/(X - ) 2 and 1/(X + ) 2 terms in (17.81) come from the two-point functions of j ± : 



(OlT^pOj^lO) 



1 + a 



(7.26) 



The resulting contribution to the cross-section corresponds to producing the stuff that is 
created out of the vacuum by the action of j ' . 



{j± stuff} 



(7.27) 



This stuff looks massless because has an integer dimension. Using the notation (I6.26p . 
the expansion (17.81) can be written as 

(0\TO(x 1 )O(x 2 )O*(y 1 )O*(y 2 )\0) 



*A W (G)*A W (C 2 ) 



i - vr 2 crc 2 



<rir> + 



Ci(3 + \d~J + )+C 2 (d-J + \j + ) , CfC 2 (d ./• d ./•} 



2 4 

cr 2 (j + |9 2 j + ) + c 2 - 2 (5 2 j+ij+) , cr 3 0' + i^j + ) + c 2 - 3 (^ J +| J +) 



6 



24 
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(i 2 &(d-j + \d*j+) + crcf(d 2 j + \d-j + ) 



12 



+ 



'1 + a) V 



Ci C2 



(T ++ |T++) 



Cr(T++|9_T ++ ) + C 2 "(9_T + +|T++) 



+ 7T C1C2 



cr(j 2 |9_j 2 ) + c 2 -(^j 2 |j 2 ) 



+ {+ -H- — for all asymmetric terms} + . . . 
This is in agreement with the OPE 
T0*(C)O(O) 



(7.28) 



6 



24 



7T 



( (c + ) 2 t- 



(r) 2 T H 



(C + ) 3 d+T- + (C) 3 d_T- 



+ n\+C(3 2 + C+d+j2+ S d - 32 + 



+ 



+ ... 



(7.29) 



where the operators on the r.h.s. are evaluated at position 0. We can verify that the contri- 
bution of the derivatives is exactly the way it must be to respect the conformal symmetry: 
the contribution of j 2 (a scalar with dimension A = 2) follows (16.30!) . An analogous ex- 
pression for the conserved current j 11 (which always has dimension D — 1 in D spacetime 
dimensions) is [9TJ [59] 

•1 



Cm 



« 



T(D/2) (-C 2 □) 



dt[t(i-t)] D/2+m -V^ 



oc 



(C 



2\d-D/2+l 



D + 2 
2 ^8(D + 1)^ ' 



( 2 B 



D + 4 . (C . 9) 3_C 2 (C-5)D 



8(D + 1) 
+ . . 



48(£> + l) v * ~' 16(D + 1) ' '"J ( ^ 7 ' 30 ' ) 

which agrees with the j M terms in (17.291) . if we take into account that in 2D the conservation 
of the vector and axial currents implies d + j + = <9_j~ = 0@ Similarly, for the stress-energy 
tensor T^ u (which always has dimension D) we have the expression [59J 



(C 



r(D/2 + i) (-C 2 a) m 

2V i-d 2+1" 4 m m\T(D/2 + 1 + m) 

m=0 ' 



dt [t(l - *)] 



D/2+m t£.0 



OC 



(C 



2 xd-£»/2+l 



1 C-9 

1 + ^— + ... 



(7.31) 



28 More precisely, in order to account for the relative signs between the j + and j terms in (|7.29[) . we 
should say that (|7.30l) describes the contribution of the axial current j ^ = —^ iV j v rather than j M . 
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which is consistent with (17.29p . 
We can also use (I6.32|) with 



X12 = C 



i • 



2^34 — C2 , ^24 — —X 



(7.32) 



to verify the contributions of the various operators directly to (I7.8p . The contribution of 
[I = 1, A = 1) should be proportional to 



iA w (d) zA w (C 2 ) [ln(l - 77-) + ln(l - V 



zA w (Ci)zA w (C 2 )__. 2 

K A ) 



-, 1 Cl ~ C2 1 Cl 2 Cl C 2 + C 2 

+ x- + (X-) 2 

(C1-C2) (cr 2 - crc 2 - + c 2 



(7.33) 



{-^+} 



which exactly reproduces the terms that are attributed to j + or j (and their derivatives) 
in (I7.28p . The contribution of T^ v (1 — 2, A = 2) is proportional to 



-12zA w (Ci)*A w (C 2 ) 
= zA w (Ci)*A w (C 2 ) 



l+(i--i)ln(l- ?7 - 



2 ._2 



Ci C2 



(x-y 



\ 2(cr - q) 

x- 



+ {-^+} 



(7.34) 



which agrees with what comes from the T ±=t terms in (17.281) . The contribution of j 2 (£ = 0, 
A = 2) is proportional to 



*A W (G)*A W (C 2 ) ln(l-7T)ln(l-77 + ) 



2/-2 



zA w (Ci)^A w (C 2 ) 



etc. 

(x 2 ) 2 



l _i_ Ci C 2 Ci" C: 



+ . . . 



(7.35) 



which agrees with the j 2 terms from (17.281) . To summarize, in terms of i"Q,A from (I6.32p . 
(17. 8p is exactly reproduced by 



l + (l + a)K 1>1 + ^^K 2 , 2 + ^±^K , 2 



(7.36) 



up to order ((/X) 5 . In order to account for even higher order terms that are not written 
explicitly in (17. 8p more operators are required. For example, we find that in order to reach 
terms of order {C/X) 11 we need to take 



1 + (1 + a) Ki,i H ^ -«- 2 , 2 H - iv , 2 H - K 3<3 H - K h3 



+ 
+ 



(1 + a) 4 (1 + a) 2 



4! 



5! 



3! 

K 1 ( 1 + a )V 1 ( 1 + a ) 4 j; 
-<m,4 i ^ -ft- 2,4 i ^ -ft-0,4 



(1 + a) 5 (1 + a) 3 



5! 



7-4! 



^5.5 + 



6 

(1 + a) 5 (1 + a) 3 



(7.37) 



4! 



5! 
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In the free-fermion limit discussed at the end of subsection 17.21 operators that have, 
besides the iA u ((), both ( + and (~ in their OPE coefficient (e.g., j 2 ), will not contribute to 
the cross-section because (17.91) vanishes for a — > unless m = or n = 0. On the other hand, 
all the operators with A = i (including j^, T M ", and others) and their derivatives contribute 
to the production of the fermions in ( 17.22p . In the free theory we have j + {x) = 2 ifjl(x)i/ji(x) , 
j~(x) = 2 ip2{ x )' l l ! 2{x) , so it is clear why the operators produce the fermion pair. However, 
these operators produce the two fermions at the same spacetime point. Processes in which 
the two internal vertices in figure fl~6k are at a finite distance from each other are described 
by the additional operators that include both derivatives of j M as well as other primary 
operators that appear in the OPE of two ip s. 



7.4 Exclusive treatment: amputated 3-point functions 



Like in section [6~5| we can verify that our inclusive results agree with exclusive computations 
that are based on amputated 3-point functions. For example, for the process 

(f) + (j) -» (f) + (j) + {j 2 stuff} (7.38) 

we know that 

(0|Tf (x)f(0)\0) = (7.39) 

«,|T/M 0> W O { y)\0) = ^f-^Zl^M^V) (7 . 40) 

A = 2, C 2 =<i±^, C 3= (1 + a »! (7.41) 



l + d-a)(|^ + |-)+... 



and then ( I6.42[) gives 
\m (i + ^) 2 [r(i-q)] 2 Q 2 . r^ 0^. r^ 2^ 

(7.42) 

in agreement with the j 2 contribution to (I7.19p . 
For the process 

<p + <p ->• <p + <p + {j» stuff} (7.43) 

we need a generalization of the analysis in section [2]to the case where the amputated operator 
is a vector. The 3-point function, that is the analog of (12. 9p . is 



-x 2 2 + ze) d - (A - 1)/2 (-x 2 3 + z e ) (A - 1)/2 (-4, + *e) (A " 1)/2 

(7.44) 



44 



The appropriate analog of the D'EPP formula (12.71) says that for Si + 62 + £3 + 1 = D [52J 

^ + 1)^2 + 1)^3 + 1) 



d D x 4 [6Z-2 



3^^3X43 v 



x 



13 



X 



14 X 24 



X 



14 X 2i 



UJ \- L 24J l x 43J 



7T 



Jb 1 Q X • 



D/2 ( x 13 _ ^23 
x 13 x 23 



+ r (f r (f - (Ji.) r (f - tfj) 



2 nD/2-5 3 / 2 \D/2-S 2 -l /■ 2 \^/2-*i-l 



12; 



(4.) 



(7.45) 



23 ' 



For solving A2.2H we use this formula as 



1 



2 sd-(A-l)/2 



X 



12/ 



rf D x 4 1^-2- 



J '43 J '43y \ / • x '14 • x '24 



-13 



A+l 



r(A + r 



a4 *24/ (x? 4 ) (A - 1)/2 (^ 4 ) (A - 1)/2 (^ 3 ) A 



7T 



D/2 / X 1Z _ 



^13 a '23/ (X 



(A-i)r(f-A) [r(*±i)] 2 

2 vd-(A-l)/2/ 2 x(A-l)/2/ 2 \(A-l)/2 



(7.46) 



12; 



[x 



13 J 



X 



23) 



Since in our case A = A, the amputated 3-point function is proportional to the ordinary 
3-point function, with the prefactor fixed by flA.14j) : 



Q ft (x 3 \x 1 ,x 2 ) 



2tt 



1 + a 



(0\TO*(x 1 )O(x 2 )f(x 3 )\0) 



The ordinary 3-point function ( 1A.16I) in momentum space is 



(<D*(-P) 0{P - Q) 3 ± {Q)) = ±2(1 + a) 



iQ~- 



Q 2 + it 

Thus the amplitude for the production of stuff is 

2tt 

Q 



[iA u {P)-iAu{P-Q)} 



Mj± = Ti^± [iA u {P) - tA u (P - Q)] 



(7.47) 



(7.4f 



(7.49) 



Alternatively, instead of using the D'EPP formula, we could do the amputation of j M 
directly in momentum space. We obtain the amputated 3-point function iQ± by dividing 
the ordinary 3-point function (17.481) by the propagator 



(r(-Q)rm = i 



,1 + a (Q±) 2 



7T Q 2 + ie 



(7.50) 



where we used ()A.14j) and (17.151) with M — 2. This again gives (I7.49p . Yet another method 
for obtaining the same result is described in appendix ID1 
From (IA.14I) and (17. 16ft the phase space of is 



Thus 



$ i± (Q) = 2(1 + a) (Q±) 2 0(Q±) 5(Q 2 ) 

")±\ a l 2 



\M J ±\ 2 $ J ± = 8n 2 (l + a) [A U (P)Y 



1-1 



(7.51) 
(7.52) 
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To compare (I7.52p with our previous results, we expand it in Q/P: 



\M j ±\ 2 $ 1 ±=8n 2 (l + a) [A u {P)f 



Q_ 

p± 

hfi 



±\ 2 



:i-a) 



Q_ 

p± 



±\ 3 



; 12 V P ± J 



(7.53) 
9(Q + )5(Q 2 ) 



This is exactly what we obtain if we repeat the computation that led to (I7.19P but include 
only terms that are attributed to ^ and its derivatives in (I7.28p . 



8 Comments about the massive bosons 

When we go to energies above the conformal region, the massive boson propagators in 
iAo(x) [see (I4.17P ]. and consequently in the higher n-point functions (I4.20p . start playing a 
role. As we have shown in [53] and reviewed in section [5j the resulting total cross-section 
for the disappearance process can be interpreted as involving phase space integrals over 
the unparticle stuff with an arbitrary number of additional massive bosons, as shown in 
figure These massive bosons have the same origin as the massive boson of the Schwinger 
model [HU |82] and are analogous to the hadrons of QCD. 

This interpretation of the disappearance process motivates us to describe the higher order 
processes of sections [6] and [7] analogously: production of the various types of unparticle stuff 
like in the conformal limit and some number of additional massive bosons. In the following 
we show that this interpretation indeed fits into a consistent picture. 



8.1 Processes with a single massive boson 




Figure 17: A term with one massive boson (the short-dashed line) in the cross-section for the missing 
energy process. 

Consider, for example, a term in which one massive boson line is attached to a dotted 
line of the missing energy process of figure [121 as shown in figure [T71 We would expect this 
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term to describe processes in which a massive boson B is produced along with unparticle 
stuff of type k: 

<P + <P^<P + <P + {O k stuff} + B (8.1) 

Naively this looks problematic because the contribution of figure [IT] to the cross-section is 
negative since similarly to (14.171) we can write the full dotted line propagator as 

1 1 1 °° (A ) n 

iA a (x) = - . . = - . . exp [ATTiaA(x)] = - . . V 1 V l iA ( x )T ( 8 - 2 ) 
uv ; iA Q (x) iA u (x) Fl v n %A u {x)^ n! 1 v n K J 

Since a < 0, the prefactors of terms that include an odd number n of massive bosons are 
negative, unlike in (I4.17P where all the prefactors are positive. However, there actually 
exist several diagrams that contribute to the process (18.11) . as shown in figure [181 including 
diagrams in which the boson is attached to a solid line. Only the sum of the four diagrams 
should be non-negative, and this is indeed satisfied since the sum is proportional tc0 

\Qk(P) - Qk(P - i)\ 2 (8.3) 

Here we use Quip) to denote the conformal 3-point function of O*, O and Ok (with the Ok 
leg amputated) in momentum space, with momentum p entering the 3-point function at the 
O* point, momentum Q leaving with the Ok stuff (it is assumed to be the same in all the 
diagrams and therefore not specified explicitly as an argument of Qk), and momentum p — Q 
leaving at the O point. 




Figure 18: The leading diagrams contributing to the cross-section of the process (18. ip . Here Qk{p) 
denotes the unparticle amputated 3-point function with incoming momentum p at the O* vertex. The 
momentum flowing through the boson line is denoted by I. The factor — Aira that is common to all the 
diagrams, as well as the common factors of h and the phase space factors, are not written explicitly. 



29 More precisely, besides the diagrams shown in figure 1181 there are other diagrams (of higher order in 
a) that contribute to the process (|8.ip . in which additional boson lines are attached between the bottom O 
and O* and/or between the top O and O* of the diagrams in figure [TH] They can be described by including 
the corresponding terms in (|8.3|i . The counting obviously works because it is the same as it would be with 
regular Feynman diagrams. 
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Figure 19: A vertex described by (j8.5|) . The long-dashed lines are the standard model particles 0, the 
solid line is the unparticle leg, and the short-dashed line is the massive boson B. The arrows indicate 
the direction of flow of the chiral charge. 



In the more common particle-physics language we can describe these processes as follows. 
We already know that in the unparticle limit (13. 19[) gives the effective interaction 

C M = \{Pu<P 2 + O u <f) (8.4) 

where Ou is the operator O in the unparticle limit (which was often denoted simply by O 
when we restricted ourselves to the unparticle limit in the previous sections). Similarly, as 
we will now verify, processes in which a single additional boson is produced can be described 
by including also the effective interaction 

Ant, i = -i ^ \ (B Ou 0* 2 - B O u 2 ) (8.5) 

where we used — Ana = (2e/m) 2 and assumed that the massive boson field B has the standard 
kinetic term ^d^Bd^B. This interaction gives the vertex V\ shown in figure [T9] (and its 
hermitian conjugate V* = — Vi). In the disappearance process with a single massive boson 
(the n — 1 diagram in figure IHJ), we have a V\ vertex and a V* vertex, which give the 
extra factor of (2e/m) 2 = —Ana relative to the unparticle-only disappearance process. This 
is consistent with (I4.27p . In the process (18. ip . the four diagrams of figure [18] contain the 
following vertices involving the massive bosons (the first factor is the vertex from the bottom 
part of the diagram and the second from the top): 

V*-Vi Vi-Vi V? -K Vx-V* (8.6) 

The products of the two vertices are positive for the first and fourth diagrams (which are 
precisely the diagrams where the boson is coming from the solid lines) and negative for the 
second and third diagrams (where the boson is coming from the dotted lines), which agrees 
with what we obtained above from our general approach. 

Even more simply, from the perspective of the effective interaction f)8.5p . the sum of the 
diagrams in figure [18] can be viewed as the standard calculation of f \A4\ 2 d§, where M. is 
the sum of the two interfering diagrams (figure 1201) that one can construct for the amplitude 
of the process (18. ip using the vertex from figure [19] and its conjugate. The relative minus 
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VtQkiP-l) V?Q k {P) 

Figure 20: The leading diagrams for the amplitude for producing one additional boson in the missing 
energy process. The circle marks the point at which the standard model particles inject the chiral charge 
into the diagram, while the square is the point at which the opposite process occurs. The massive boson 
lines are short-dashed and carry momentum £. 



sign in ( 18. 3ft appears because one diagram in figure EO has a factor of V\ while the other has 

V{ = -v x . 

It is interesting to note that there exists a missing energy process in which only the 
massive boson is produced, without the unparticle stuff: 

<J) + <j) ^ <j) + <J) + B (8.7) 

This arises similarly to (18.11) . except that O k is the identity operator in (17.291) which gives 
rise to the disconnected term in (17.281) and (17. 8p . This is different from the disappearance 
process of section [5] where the massive bosons can only be produced together with unparticle 
stuff. 

We can also approach the process (18.71) . and similarly the processes (18. ip . from a different 
point of view by considering the operator jg (or equivalently its axial counterpart j^ 1 = 
' 3Su) that is discussed appendix [A] In particular, the OPE contains the term 

T<T(C)O(0) D -27n*A w (C)C,J5 M (0) (8.8) 
The form of the two-point function of jg, (1A.17I) - (1A.18I) . suggests the replacement 

5 f^> fj±PB (8.9) 



7T 

where we use 

(□ +m 2 )iA(x) = -iS 2 (x) (8.10) 

to reconcile (1A.18j) (up to a contact term, which is ambiguous anyway). Then we can write 

T 0* (C) O(0) D z^ tAu(C) C^B(O) (8.11) 
This means that for small y — x 

T 0*(y) 0(x) D iA u (y - x) [B{y) - B{x)\ (8.12) 

which is precisely proportional to the amplitude we would write for the process ( 18. 7p by 
summing the two diagrams in figure | 
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Figure 21: The leading position-space diagrams for the amplitude for producing just a boson, without 
the unparticle stuff, in the missing energy process. 

8.2 Processes with multiple massive bosons 



v„ = [i-Yh 
m j 



✓ \ 



Figure 22: An n-boson vertex described by (j8. 13j) . The long-dashed lines are the standard model 
particles <f>, the solid line is the unparticle leg, and the short-dashed lines are the n massive bosons B. 
The arrows indicate the direction of flow of the chiral charge. 



The ideas discussed above can be generalized to processes with an arbitrary number of 
bosons by writing 



•-int 



O u (j)* z + e i{2e/m)B 0* u cj) 2 ) 



(8.13) 



This expression contains (18. 4p . (18. 5p . and similar multiple-boson vertices, as shown in fig- 
ure |22j One can see that (18.131) is equivalent to our original description by noticing that the 
interaction (I3.19P gives 

(0\T<p 2 (x 1 )... ( p 2 (x n ) ( p* 2 (y 1 )...^ 2 (y n )\0) 
oc (0\TO(x 1 ) . . . 0{x n ) 0\y x ) . . . O*(y n )\0) 



llj,fc ex P [(2e/m) 2 iA(xj 



Vk 



rifc>j ex P [(2e/m) iA(xj - x k )] exp [(2e/m) zA(% - y k )] 



x(0\TO u (x 1 ) . . . O u (x n ) 0* u { Vl ) . . . 0* u {y n )\ti) 
where we used ( fODjl with ( H7T7|t . while (J8TT3]) gives 



(8.14) 



(0\T ( j ) '(x 1 )... ( j ) '(x n ) ( j ) *\y 1 )...cP*\y n )\0) 
oc (0|Te~ i(2e/m)B{a;i) ...e- i 



(2e/m)B(x n ) J(2e/m)B( yi ) 



. . . e 



i(2e/m)B(y n ) 



|0> 



x (0\TO u ( Xl ) . . . O u [x n ) 0* u ( yi ) . . . O* u (y n )\0) 



(8.15) 
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After contracting the boson operators in (I8.15P we indeed obtain the result (I8.14p . 





ix'l 

1 
1 
1 


X 

X 

X 


1 
1 
1 

© 


1 
% 
X 

□ 




Q k (P 



Qk{P) 



-Qk(P 



-Qk(P 



Figure 23: The leading diagrams for the amplitude for producing 2 additional bosons in the missing 
energy process. For clarity, only the boson lines are shown explicitly (but the unparticle 3-point function 
is also assumed to be present like in figure I2TJ1) . Common prefactors are not written. 



To see more explicitly how (I8.13P works, let us look at missing energy processes that 
produce 2 bosons. Using the vertices V\ and from figure [22] (and their conjugates), the 
amplitude for such a process is given by the sum of the diagrams in figure [23] Squaring the 
amplitude and including a factor of 1/2 due to the 2 identical bosons in the final state, the 
cross-section is proportional to 

' f d$ \Q k (P - l x - i 2 ) + Q k (P) -Q k (P- £x) -Q k (P- £2) | 2 (8. 16) 



On the other hand, this is precisely the result we get from our general formalism by summing 
the diagrams in figure [22] where we used (I4.17P and (18.21) for computing the effects of the 
bosons] 30 ! 

As another check of (18. 13ft . consider the case in which n boson lines are attached between 
the two ends of an internal solid or dotted line (and for simplicity, no other bosons are 
attached to these vertices). From (I4.17P and (18.21) we have a factor of 

(-4™)™ (4vra)" 

r 1 - or ^-L 8.17 

n! n\ 

in the case of a solid or dotted line, respectively. On the other hand, taking the product of 
the two relevant vertices (using figure [22] and/or its conjugate) gives 



2e 
±i— 
m 



.2e 
m 



-Ana) 1 



or 



±i— ^ (±z-l =(4vra)" 



Hi 



for these two cases. Dividing the expressions in (18.181) by the symmetry factor n! we obtain 
an agreement with ( 18.171) . 



8.3 Massive boson decay 

At the leading order in h, the massive boson B decays as 



.19) 



30 To make the agreement manifest we need to rename the integration variables i\ f-> £2 in some of the 
terms that we obtain after expanding the square in (I8.16[) . 
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Figure 24: The leading diagrams for the cross-section of producing 2 additional bosons in the missing 
energy process. For clarity, only the boson lines are shown explicitly. Common prefactors and phase 
space factors are not written. 



via the vertex in figure [TH] and its conjugate. The rate for each of these two modes is given 
by 

r = ^l-M| 2 $ (8.20) 

where 

\M\ 2 = %h 2 (8.21) 

and 

A(a) f d 2 p u 



(£m) 2a J (2tt) 



x 



n 



d 2 Pi 



2ir5(p 2 -m 2 )6(p°) 



(2vr) 5 (p B -pu-Pi- P2) 
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A(a) ^(iV 1 -^ 2 ) dpf_ 



L 



~ pt dpi 




) 



a 



(8.22) 



At the leading order in a and for Cmwe obtain 



r ~ 




7T 



12 



(8.23) 



9 Conclusions 

The Sommerfield model is a useful toy model of the Banks-Zaks sector. It becomes scale- 
invariant in the infrared, with fractional anomalous dimensions. Since all of its correlation 
functions can be computed exactly, we can answer explicitly many of the questions regarding 
the physics of the Banks-Zaks sector as seen by the (toy) standard model observer, for an 
arbitrary coupling strength in the conformal sector. In particular, we were able to explore the 
behavior of the theory away from the low-energy scale-invariant regime and to incorporate 
the unparticle self-interactions. Most importantly, we believe, we used this toy model to 
provide consistency checks for the two new ideas in this paper: the extension of unparticle 
phase space calculations beyond the two point function to processes involving unparticle self- 
interactions; and the interpretation of the result using the conformal partial-wave expansion 
in terms of an "amplitude" for production of different types of unparticle stuff. These results 
are sufficiently general that they should apply to unparticle theories in 3+1 dimensions. 
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A Operators in the Sommerfield model 

The description of the Sommerfield model in terms of the free fields V, A in (13. 5 p imme- 
diately allows us to compute any correlation function involving these fields. Using this fact 
we were able to compute arbitrary correlation functions of the original interacting fermionic 
fields ij} in section 13.11 and it is straightforward to extend these results to correlation func- 
tions that involve A 11 as well. We have also analyzed in section H] the composite operator O 
that appears in the OPE 



TV 2 *(C)^i(0) 



Co(C) 



O(0) + --- 



(A.l) 



c(0 



53 



In this appendix we would like to discuss the properties of additional operators that exist in 
the theory, obtained using similar methods. In particular, let's consider the vector current 
that is classically given by = tfi^ip. We expect to find j + in the OPE 



T^i(C)^(0) = Te- fe Ma + ^ c »tfi(C)e*< v (°) + ^ o »*I(0) 



= ^(0 + 



cb(C)c(0 



J + (0)-^(d- + ^d + )(V(0)+A(0)) 



where J + is the free current from 



TM/ 1 (C)^(0)=^ 1 (C) + -J + (0) + 



+ ... (A.2) 



(A.3) 



The auxiliary fields V and A are our own constructions. In the theory, they appear only in 
the combination A 1 . Thus we see from here (and the analogous expression for Ti/^O^sKO)) 
that we should define 

(A.4) 



2e 



r = J T ~d T A 

71 



which gives 



T^ 1 (CM(o) = ^ 1 (C) + 



t^ 2 (C)^(o) = ^ 2 (C) + 



c (C)c(C) 



C (QC(Q 



j + (o)-^(a + (o) + ^-(o) 

r(o)-^(^ + (o) + A-( ) 



+ . . . 



+ . . . 



(A.5) 



(A.6) 



It can be shown by deriving Ward identities that both j M and A 1 are conserved vector 
currents, but they are not independent: 

(0|T^(x)A"(x)|0)^0 (A.7) 

It is useful to replace j M and A 1 by two linear combinations that do not mix. These are 

j^ = (l + a)(r--AA and f s = -af + (1 + a)-A» (A.8) 

V 71 / 7T 

Rewriting the OPEs in terms of and jg we have 



t^ 1 (C)Vi*(o) = ^ 1 (C) + 



Co(C)c(C) 



2 + a 
1 + a 



#(0) 



-a C 
1 + aC 



Jt(0) 



+ ... 



(A.9) 



T^ 2 (C)^(0)=^ 2 (C) + 



Co(C)c(C) 



2 + a 
1 + a 



Jt(0) 



-a C" 



l + aC + 



#(0) 



+j'J(o)-^^(o) 



+ ... (A.10) 
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The current jj, happens to be exactly the combination whose axial counterpart = 
~ ef11 ' 3tv is anomaly-free (unlike j 5fl = —e^jj) UM- ^ nas the non-canonical ratio of (1 + a) 
between the axial and vector charge of ipi, exactly as in the Thirring model [69j EH ED 192] - 
The current jg can be dropped at low energies because its correlation functions always involve 
the massive propagator A(x) from ( 13.141) . As a result, at low energies the OPEs reduce to 
those of the Thirring model where j£ is the only current F^l The OPE ( 1A.1I) also agrees with 
the Thirring model at low energies [70] . In terms of our parameters, the Thirring model 
interaction is 



na , 2 e 2 .,, 



£ «=2(T7^ = -^ (A ' n) 

where j\ = j^jr ^ — 3t3t- A more general treatment of the operators in the Thirring model 
can be found in [89] . 

For the product 0(x)0*(y) to which we refer in section ITUl we find the OPE 

T0*(C) O(0) = tA o (0 [1 + 2 7 ne / a / T(0) - 2tt 2 (e^ff (0) 

+*ie>» ((%U0) + (»( p d„f(0)) +...} (A.12) 

In the unparticle limit we can again drop jg to obtain 

T 0*(() O(0) = iA u (() [ 1 + «r [( + j T (0) - C#(0)] 

+ 4 [(C + ) 2 d + Jr(0) - (C) 2 d-Jt(0)] (A.13) 



7T 



2 



(C + ) 2 J^ 2 (0) + (C) 2 jf(0) + vr 2 C + C"j|(0) + • • • 



where we used the fact that since both jj, and j^ 1 are conserved d±j^(x) = 0. The two-point 
functions of are§ 



(0|Tj±(*)j±(0)|0) = - i y (0|Tj+(x)j-(0)|0) = (A.14) 



7T 2 (x 2 - ie) 



The x-dependence is actually fixed by the requirement that j£ cannot have an anomalous 
dimension since it generates a symmetry. The operators j^ 2 are the components T ±:t of the 
stress-energy tensor [71] : 



31 This is the reason for the name jt- Similarly, js gets its name because it is the only current that is 
present in the Schwinger model. 

32 It is worth saying here that (0|Tjy (x)j^ (0) 1 0) may contain a contact term which is ambiguous since 
current 2-point functions in 1 + 1 dimensions are subject to subtractive renormalization. This is analog in 
1 + 1 of the ambiguity of current 3-point functions in 3 + 1 that leads, for example, to different forms for the 
anomaly. 
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while j\ is a scalar operator. The correlation functions of such products of currents can be 
obtained (apart from the contact terms) by contracting the currents like free fields. It is also 
useful to know the 3-point function 

(0\TO*( Xl ) <D(x 2 )j±(x 3 ) |0> = T (l + a)- (-^-. iA u {x 12 ) (A.16) 

The two-point function of the current j$, which we discuss further in section 18.11 is 

(0|Tj|(x)j|(0)|0) = -dizA(x) (A.17) 



7T 

(0\Tjt{x)jg(0)\0) = -m 2 iA(x) - a(l + a) - 5 2 (x) (A.18) 

7T 7T 

B Basic properties of lightcone coordinates 

The lightcone coordinates x ± = x° ± x 1 in 2D have the properties 

x± = ^- x 2 = x + x- d± = \ (d ± dx) □ = 4 <9+<9_ (B.l) 

9+- = 9-+ = \ 9 + ~ = <T + = 2 e- + = -e + - = 2 e_ + = -e + _ = ~ (B.2) 

d 2 p = ^dp + dp- 6(p°) d(p 2 ) = d(p + ) 6(p-) (B.3) 

Under Lorentz transformation, components of vectors transform as A ± — > e TV A ± , where 
i] is the rapidity, so objects of the form A + B~, A + /B + , A~/B~ are scalars. They can be 
written in the more conventional Lorentz-invariant form as 

A+B- = A U B» + (TA U B V 4x = , „ ^ ; „ (B.4) 

C Numerical treatment of IR divergences 

In this appendix we consider the simple (unphysical) diagram in figure [25] in order to show 
explicitly how to treat the IR divergences of the dotted line in momentum space. 
In position space, the diagram is simplyr^l 

%M = - iA a {xi - x 2 ) iA a {xi - x 2 ) = -1 (C.l) 

so in momentum space M. = i(2ir) 2 5 2 (P), which gives the "phase space" 

$ = 2(2n) 2 6 (P°)5 2 (P) (C.2) 



33 Since this is not a diagram that actually exists in our theory, we do not expect factors of i to work out, 
so we added a minus sign in order for it to match with what we get using our algorithm. 
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h = p 



Figure 25: As before, the solid line represents the propagator iA and the dotted line iA. The 
momentum flowing into the diagram is P = (M, 0). 

On the other hand, if we want to do the calculation in momentum space, as we are suggesting 
to do for the diagrams of section I7.1[ we need to compute the integrao 

/= J d%<f£ 2 9(f l )9(£l)F(e 1 )9(el)B(£l)(£l) a S 2 (P-e 1 -£ 2 ) (C.3) 

According to one method of calculation, we take 

F{h) = h(£i) (C.4) 

This also needs to differentiated with 

d 2 

(C.5) 



in order to properly account for the dotted line. But using the 5-function to write l\ = P — £2 
we can just differentiate with 

<9 2 8 B 2 

P 2 J^ (C6) 



dP+dP- d(P 2 ) d(P 2 ) 2 
after computing the integral jffl According to the second method we just take 

m) = fM+%fx(fi) (c.7) 

The resulting $(P 2 ) should reproduce the 5-function from flC.2j) in the sense that 



dP + / dP~ -> 1 (C.8) 

Jo 2(2vr) 2 1 ; 

for Pmax ^> A. The integral is finite despite the fact that $(P 2 ) is a function of P 2 alone 
because 

p2 



max 

>2 /R/ 7)2 



dP 2 $(P 2 )^0 (C.9) 







34 0ne should further multiply this by 16 sin 2 (to) in order to obtain the phase space with the correct 
prefactor. 

35 We used the regulated f\(l\) in (|C4[) for numerical convenience. 
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Figure 27: The solid line presents the two-dimensional integral over our numerical "<5-function" , as in 
(|C.8ft . The dashed line (normalized arbitrarily) demonstrates (jC.9j) . (For -P^ax* we use un ' ts ' n which 
A 1.) 



To confirm the viability of this procedure, we have performed the calculations described 
here numerically [90]. We used 

with a regulated ^-function 

3 

^rcg(^) = c n exp(-nx) (C.ll) 

n=0 

with the coefficients c n chosen such that Teg (0) = 9' Te JO) = and reg (x) — >■ 1 for x ^> 1. 
The results for a = —1/2 are presented in figures 1261 and [271 The specific data we show were 
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obtained using ( 1C.7j) . A computation using (1C.4|) yielded results that were essentially the 
same. 

D An alternative method for amputating the leg 

Based on the observation that for j M we have A = A, we expect (and confirm in the following) 
that the amputated 3-point function is proportional to the ordinary 3-point function. The 
latter is given by (IA.16j) and can be written as 

(0\Tf(x)O*{x 1 )O(x 2 )\0) = e^dJ(x\x 1 ,x 2 ) (D.l) 

where the derivative is with respect to x and 

f(x\x 1 ,x 2 ) = (1 + a)-^- [ln(-0 - x 2 ) 2 + ie) - ln(-(x - Xl ) 2 + ie)] iA u (x 1 - x 2 ) (D.2) 



We also write (IA.14D as 



(0\Tf(x)f(0)\0) = j±l (d^\n(-x 2 + ie) - 2mg^6 2 (x)) (D.3) 



For using ( \2.2\i , note that 



J" = / d 2 xd^d v \n(-(x' -x) 2 + te)e up d p \n(-(x-x 1 ) 2 + ie) 



\2 



= Ame^d v \n(-(x' - x^ 2 + te) (D.4) 
Then the effect of the first term of flD.3j) is to multiply the 3-point function by 

(1 + a)- (D.5) 

7T 

At the same time, the second term of (ID. 31) multiplies the 3-point function by 

-(l + o)^- (D.6) 
Thus the amputated 3-point function that would satisfy 

d 2 x (0|TyV):T(x)|0) iQ„{x\xi,x 2 ) = (0|T^(x')C*(xi)C(x 2 )|0) (D.7) 



is 

Q»(x\x u x 2 ) = (OlTjVWixJOfaW) (D.8) 

1 + a 

which agrees with (17.47j) . 



59 



References 



[1] H. Georgi, "Unparticle physics," \Phys. Rev. Lett. 98 (2007) 221601 
|arXiv : hep-ph/0703260| 

[2] H. Georgi, "Another odd thing about unparticle physics," 



\Phys. Lett. B650 (2007) 275HarXiv: 0704. 2457 [hep-ph] [ 

[3] A. A. Belavin and A. A. Migdal, "Calculation of anomalous dimensionalities in 
non-Abelian field gauge theories," Pisma Zh. Eksp. Teor. Fiz. 19 (1974) 317. 



JETP Lett. 19 (1974) 181 



[4] W. E. Caswell, "Asymptotic behavior of non-Abelian gauge theories to two-loop 
order," \Phys. Rev. Lett. 33 (1974) 244[ 

[5] T. Banks and A. Zaks, "On the phase structure of vector-like gauge theories with 



massless fermions," Nucl. Phys. B196 (1982) 189. 



M. J. Strassler, "Why unparticle models with mass gaps are examples of hidden 
valleys," |arXiv: 0801 .0629 [hep-ph] | 

[7] B. Grinstein, K. Intriligator, and I. Z. Rothstein, "Comments on unparticles," 



Phys. Lett. B662 (2 008) 367] |arXiv: 0801 . 1140 [hep-ph] 



[8] N. Seiberg, "Electric-magnetic duality in supersymmetric non-Abelian gauge 



theories," Nucl. Phys. B435 (1995) 129, arXiv:hep-th/9411149 



[9] P. J. Fox, A. Rajaraman, and Y. Shirman, "Bounds on unparticles from the Higgs 
sector," \Phys. Rev. D76 (2007) 075004| |arXiv: 0705. 3092 [hep" = ph]] 

[10] Y. Nakayama, "Supersymmetry breaking effects in the unparticle sector and conformal 



sequestering," Phys. Rev. D76 (2007) 105009, arXiv: 0707. 2451 [hep-ph] 



[11] T. A. Ryttov and F. Sannino, "Conformal windows of SU(iV) gauge theories, higher 
dimensional representations and the size of the unparticle world," 
\Phys. Rev. D76 (2007) 105004[ |arXiv : 0707. 3166 [hep-th] 



[12] T. A. Ryttov and F. Sannino, "Supersymmetry inspired QCD beta function," 
\Phys. Rev. D78 (2008) 0650011 |arXiv: 071 1.3745 [hep-th] | 

[13] F. Sannino, "Dynamical stabilization of the Fermi scale: Phase diagram of strongly 
coupled theories for (Minimal) Walking Technicolor and unparticles," 



arXiv: 0804. 0182 [hep-ph] 



[14] J. Polchinski, "Scale and conformal invariance in quantum field theory," 
\Nucl. Phys. B303 (1988) 226| 



60 



[15] D. Dorigoni and S. Rychkov, "Scale invariance + unitarity =>- conformal invariance?," 
|arXiv: 0910 .1087 [hep-th] | 

[16] G. Mack, "All unitary ray representations of the conformal group SU(2,2) with 



positive energy," Commun. Math. Phys. 55 (1977) 1. 



[17] I. T. Todorov, M. C. Mintchev, and V. B. Petkova, Conformal invariance in quantum 
field theory. Scuola Normale Superiore, Pisa, 1978. 

[18] D. M. Hofman and J. Maldacena, "Conformal collider physics: Energy and charge 
correlations," JEEP 05 (2008) 012[ |arXiv : 0803 . 1467 [hep-th]] 

[19] M. Bander, J. L. Feng, A. Rajaraman, and Y. Shirman, "Unparticles: Scales and high 



energy probes," \Phys. Rev. D76 (2007) 11 5002, a rXiv: 0706 .2677 [hep-ph] 



[20] T. Kikuchi and N. Okada, "Unparticle physics and Higgs phenomenology," 
\Phys. Lett. B661 (2008) 360HarXiv: 0707 .0893 [hep-ph] [ 

[21] A. Delgado, J. R. Espinosa, and M. Quiros, "Unparticles-Higgs interplay," 
\JEEP 10 (2007) 0941 |arXiv : 0707 . 4309 [hep-ph] [ 

[22] J. M. Maldacena, "The large N limit of super conformal field theories and 



supergravity," Adv. Theor. Math. Phys. 2 (1998) 231, arXiv:hep-th/9711200, 



[23] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, "Gauge theory correlators from 
non-critical string theory," \Phys. Lett. B428 (1998) 105[ |arXiv : hep-th/9802109] 

[24] E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2 (1998) 
253, [arXiv : hep-th/9802150| 

[25] L. Randall and R. Sundrum, "An alternative to compactification," 
\Phys. Rev. Lett. 83 (1999) 4690] |arXiv : hep-th/9906064] 

[26] L. Randall and R. Sundrum, "Large mass hierarchy from a small extra dimension," 



Phy s. Rev. Lett. 83 (19 99) 3370] |arXiv : hep-ph/990522l] 



[27] K. Cheung, W.-Y. Keung, and T.-C. Yuan, "Collider signals of unparticle physics," 



Phys. Rev. Lett. 99 (2007) 051803, [arXiv: 0704. 2588 [hep- ph] , 



[28] Y. Liao, "Bounds on unparticle couplings to electrons: from electron g — 2 to 



positronium decays," \ Phys. Rev. D76 (20 07) 056006] |arXiv: 0705 .0837 [h ep-ph] 



[29] S. Zhou, "Neutrino decays and neutrino-electron elastic scattering in unparticle 



physics," \Phys. Lett. B659 (20 08) 336] |arXiv: 0706 .0302 [hep-ph] 



[30] G. Cacciapaglia, G. Marandella, and J. Terning, "The AdS/CFT/unparticle 
correspondence," | JEEP 02 (2009) 049[ |arXiv : 0804 . 0424 [hep" z ph]] 



61 



[31] 

[32] 

[33] 

[34] 

[35] 

[36] 

[37] 

[38] 
[39] 

[40] 

[41] 

[42] 

[43] 

[44] 

[45] 

[46] 



M. Perez- Victoria, "Randall-Sundrum models and the regularized AdS/CFT 



correspondence," JEEP 05 (2001) 064, arXiv:hep-th/0 105048 



M. Perez- Victoria, "Unparticle signals with a few particles," \JEEP 01 (2009) 011[ 
|arXiv : 0808 . 4075 [hep-ph] | 

C. M. Ho and Y. Nakayama, "Unparticles and holographic renormalization group," 
\JEEP05 (2009) 0811|arXiv: 0903. 0420 [hep-th]] 

A. Friedland, M. Giannotti, and M. Graesser, "On the RS2 realization of unparticles," 



Phys . Lett. B678 (2009) 149[ |arXiv: 0902 . 3676 [hep-th] 



G. Cacciapaglia, G. Marandella, and J. Terning, "Colored unparticles," 
\JEEP 01 (2008) 070[ |arXiv : 0708 . 0005 [hep-phj 



Y. Liao, "Some issues in a gauge model of unparticles," 

\Eur. Phys. J. C60 (2009) 125-134HarXiv: 0804 .4033 [hep-ph] | 

J. Galloway, J. McRaven, and J. Terning, "Anomalies, unparticles, and Seiberg 
Duality," \Phys. Rev. D80 (2009) 105017[ |arXiv : 0805 . 0799 [hep" = ph]] 



A. L. Licht, "The gauged unparticle action," arXiv: 0806 .3596 [hep-th] 



A. Ilderton, "Unparticle actions and gauge invariance, 



Phys. Rev. D79 (2009) 025014, arXiv: 0810. 3916 [hep-th] 



M. A. Stephanov, "Deconstruction of unparticles," Phys. Rev. D76 (2007) 035008 
|arXiv : 0705 . 3049 [hep-ph] | 



D. Stancato and J. Terning, "The Unhiggs," JEEP 11 (2009) 101 



arXiv: 0807. 3961 [hep-ph] 



A. Falkowski and M. Perez- Victoria, "Holographic Unhiggs," 



Phy s. Rev. D79 (2009) 035005[ |arXiv: 0810 .4940 [hep-ph]] 



P. Kovtun and A. Ritz, "Black holes and universality classes of critical points," 
\Phys. Rev. Lett. 100 (2008) 171606[|arXiv: 0801 .2785 [hep-th] [ 

J. Erlich, E. Katz, D. T. Son, and M. A. Stephanov, "QCD and a holographic model 



of hadrons," \Phys. Rev. Lett. 95 (2005) 261602[ |arXiv : hep-ph7b 501 128 , 

L. Da Rold and A. Pomarol, "Chiral symmetry breaking from five-dimensional 



spaces," \Nucl. Phys. B721 (2005)T9] |arXiv:hep-ph/0501218] 

S. S. Gubser, "AdS/CFT and gravity," \Phys. Rev. D63 (2001) 084017] 



arXiv : hep-th/9912001 



62 



[47] 

[48] 

[49] 
[50] 

[51] 

[52] 

[53] 

[54] 

[55] 

[56] 

[57] 
[58] 

[59] 

[60] 
[61] 



N. Arkani-Hamed, M. Porrati, and L. Randall, "Holography and phenomenology," 



J HEP 08 (2001 ) 017] |arXiv : hep-th/0012148 



J. L. Feng, A. Rajaraman, and H. Tu, "Unparticle self-interactions and their collider 



implications," Ph ys. Rev. D77 (2008) 075007[ |arXiv: 0801. 1534 [hep-ph] | 



P. H. Ginsparg, "Applied conformal field theory," arXiv:hep-th/9 108028 



H. Osborn and A. C. Petkou, "Implications of conformal invariance in field theories for 



general dimensions," \Ann. Phys. 2 31 (1994) 311, arXiv : hep-th/930 7010 



E. S. Fradkin and M. Ya. Palchik, "Recent developments in conformal invariant 
quantum field theory," \Phys. Rept. 44 (1978) 249 



E. S. Fradkin and M. Ya. Palchik, "New developments in D-dimensional conformal 
quantum field theory," \Phys. Rept. 300 (1998) Tj 

H. Georgi and Y. Kats, "Unparticle example in 2D," 



Phys . Rev. Lett. 101 (2008) 131603HarXiv: 805 .3953 [hep-ph] 



S. Ferrara and G. Parisi, "Conformal covariant correlation functions," 
\Nucl. Phys. B42 (1972) 28lj 

S. Ferrara, A. F. Grillo, G. Parisi, and R. Gatto, "Covariant expansion of the 



conformal four-point function," Nucl. Phys. B49 (1972) 77 



V. K. Dobrev, V. B. Petkova, S. G. Petrova, and I. T. Todorov, "Dynamical derivation 
of vacuum operator product expansion in Euclidean conformal quantum field theory," 
\Phys. Rev. D13 (1976) 8871 

G. Mack, "Convergence of operator product expansions on the vacuum in conformal 



invariant quantum field theory," Commun. Math. Phys. 53 (1977 ) 155| 



E. D'Hoker, S. D. Mathur, A. Matusis, and L. Rastelli, "The operator product 
expansion of = 4 SYM and the 4-point functions of supergravity," 
\Nucl. Phys. B589 (2000) 38HarXiv:hep-th/99 11222 



G. Arutyunov, S. Frolov, and A. C. Petkou, "Operator product expansion of the 
lowest weight CPOs in M = 4 SYM 4 at strong coupling," 



Nucl. Phys. B586 (2000) 547, arXiv:hep-th/0005182 



F. A. Dolan and H. Osborn, "Conformal four point functions and the operator 
product expansion," \Nucl. Phys. B599 (2001) 459[ [arXTv:hep-th/00Tl040, 



M. D'Eramo, L. Peliti, and G. Parisi, "Theoretical predictions for critical exponents at 
the A-point of Bose liquids," \Lett. Nuovo Cim. 2 (1971) 878| 



63 



[62] 

[63] 

[64] 

[65] 

[66] 

[67] 

[68] 
[69] 
[70] 

[71] 

[72] 

[73] 

[74] 

[75] 

[76] 

[77] 



D. E. Diaz and H. Dorn, "On the AdS higher spin / O(N) vector model 



correspondence: Degeneracy of the holographic image," JHEP 07 (2006) 022 
arXiv : hep-th/0603084| 



C. M. Sommerfield, "On the definition of currents and the action principle in field 
theories of one spatial dimension," Ann. Phys. 26 (1964) 1| 

L. S. Brown, "Gauge invariance and mass in a two-dimensional model," 
\Nuovo Cim. 29 (1963) 6171 

W. E. Thirring and J. E. Wess, "Solution of a field theoretical model in one space- one 



time dimension," Ann. Phys. 27 (1964) 331, 



D. A. Dubin and J. Tarski, "Interactions of massless spinors in two dimensions," 



\Ann. Phys. 43 (1967) 263 



C. R. Hagen, "Current definition and mass renormalization in a model field theory,' 
Nuovo Cim. 51A (1967) 1033} 



J. S. Schwinger, "Gauge invariance and mass. II," Phys. Rev. 128 (1962) 2425 



W. E. Thirring, "A soluble relativistic field theory," Ann. Phys. 3 (1958) 91 



K. G. Wilson, "Operator-product expansions and anomalous dimensions in the 



Thirring model," Phys. Rev. D2 (1970) 1473 



H. Georgi, "Energy-momentum tensors and scale invariance in the Thirring model," 
\Phys. Rev. D2 (1970) 2908] 

R. Roskies and F. Schaposnik, "Comment on Fujikawa's analysis applied to the 



Schwinger model," Phys. Rev. D23 (1981) 558 



H. Georgi and J. M. Rawls, "Anomalies of the axial-vector current in two dimensions," 
\Phys. Rev. D3 (1971) 874 



J. H. Lowenstein and J. A. Swieca, "Quantum electrodynamics in two- dimensions," 
\Ann. Phys. 68 (1971) 172] 

G. Segre and W. I. Weisberger, "Investigations in two-dimensional vector-meson field 
theories," \Phys. Rev. DIP (1974) 1767} 

C. R. Hagen, "Unification of the soluble two-dimensional vector coupling models," 



Phys. Rev. D55 (1997) 10~2Tj |arXiv : hep-th7 9409038 , 



C. G. Callan, R. F. Dashen, and D. J. Gross, "The structure of the gauge theory 



vacuum," Phys. Lett. B63 (1976) 334 



64 



79 



SO 



82 

83 
84 

85 
86 



89 

90 
91 



S. R. Coleman, R. Jackiw, and L. Susskind, "Charge shielding and quark confinement 



in the massive Schwinger model," \Ann. Phys. 93 (1975) 267 



S. R. Coleman, "More about the massive Schwinger model,' 



Ann. Phys. 101 (1976) 239, 



S. R. Coleman, "There are no Goldstone bosons in two dimensions," 

3 



\Commun. Math. Phys. 31 (1973) 259, 



K. G. Wilson, "Non-Lagrangian models of current algebra," 
\Phys. Rev. 179 (1969) 14991 

A. Casher, J. B. Kogut, and L. Susskind, "Vacuum polarization and the absence of 



free quarks," Phys. Rev. DIP (1974) 732 



I. M. Gel'fand and G. E. Shilov, Generalized functions, vol. 1. Academic Press, 1964. 

A. Casher, J. B. Kogut, and L. Susskind, "Vacuum polarization and the quark-parton 
puzzle," \Phys. Rev. Lett. 31 (1973) 792| 

A. Petkou, "Conserved currents, consistency relations, and operator product 
expansions in the conformally invariant O(N) vector model," 



Annals Phys. 249 (1996) 180, arXiv:hep-th/9410093. 



S. Ferrara, A. F. Grillo, and R. Gatto, "Manifestly conformal covariant 



operator-product expansion, " \Lett. Nuovo Cim. 2 (197ijl363 



S. Ferrara, A. F. Grillo, and R. Gatto, "Tensor representations of conformal algebra 
and conformally covariant operator product expansion," Ann. Phys. 76 (1973) 161] 

R. Rattazzi, V. S. Rychkov, E. Tonni, and A. Vichi, "Bounding scalar operator 
dimensions in 4D CFT," | JEEP 12 (2008) 031[ |arXiv : 0807 . 0004 [hep" z th]] 

M. Liischer, "Operator product expansions on the vacuum in conformal quantum field 
theory in two space-time dimensions," Commun. Math. Phys. 50 (1976) 2~3~] 

Wolfram Research, Inc., "Mathematica, Version 6.0." Champaign, IL, 2007. 

K. Lang and W. Riihl, "The critical 0(N) a- model at dimensions 2 < al < 4: a list of 
quasiprimary fields," \Nucl. Phys. B402 (1993) 573| 



[92] G. F. Dell' Antonio, Y. Frishman, and D. Zwanziger, "Thirring model in terms of 
currents: Solution and light cone expansions," Phys. Rev. D6 (1972) 988, 



65 



